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Abstract: Current astrophysical research suggests that the most persistently luminous objects 
in the Universe are powered by the flow of matter through accretion disks onto black holes. 
Accretion disk systems are observed to emit copious radiation across the electromagnetic 
spectrum, each energy band providing access to rather distinct regimes of physical conditions 
and geometric scale. X-ray emission probes the innermost regions of the accretion disk, 
where relativistic effects prevail. While this has been known for decades, it also has been 
acknowledged that inferring physical conditions in the relativistic regime from the behavior 
of the X-ray continuum is problematic and not satisfactorily constraining. With the discovery 
in the 1990s of iron X-ray lines bearing signatures of relativistic distortion came the hope 
that such emission would more firmly constrain models of disk accretion near black holes, 
as well as provide observational criteria by which to test general relativity in the strong 
field limit. Here we provide an introduction to this phenomenon. While the presentation is 
intended to be primarily tutorial in nature, we aim also to acquaint the reader with trends in 
current research. To achieve these ends, we present the basic applications of general relativity 
that pertain to X-ray spectroscopic observations of black hole accretion disk systems, 
focusing on the Schwarzschild and Kerr solutions to the Einstein field equations. To this we 
add treatments of the fundamental concepts associated with the theoretical and modeling 
aspects of accretion disks, as well as relevant topics from observational and theoretical X-ray 
spectroscopy. 

PACS Nos.: 32.30.Rj, 32.80.Hd, 95.30.Dr, 95.30.Sf, 95.85.Nv, 97.10.Gz. 97.80.Jp, 98.35.Mp, 
98.62.Mw 



1. Introduction 

The publication of Einstein's General Theory of Relativity (GR) in 1915 was the culmination of 
several years of effort that began as early as 1907, when Einstein began to address the incompatibility 
between the Special Theory of Relativity and Newtonian gravitation. The nature of the inconsistency 
is clear For example, Newton's formulation of gravitation is embodied in the Poisson equation for the 
gravitational potential, W'^cf) — AirGp, where the source term p is the mass density, and it is implied 
that the field responds instantaneously to changes in p, thereby violating special relativity. To Einstein, 
this inconsistency foreshadowed a complete revision of gravitation theory. 

For a relativistic theory of gravity, the characteristic length scale for physical effects near a mass 
M, on dimensional grounds, is GM/c?, which is known as the gravitational radius. In assessing the 
magnitude of GR effects in the vicinity of a mass M, at a distance r from the mass, one compares 
r with this characteristic length scale. Therefore, we expect only small perturbations to Newtonian 
physics whenever GM /rc^ <^ 1. This is just the case for the first proposed tests of GR, and in fact, for 
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all possible terrestrial and solar system experiments. Most notable were the problem of the precession 
of the perihelion of Mercury and the prediction of the bending of starlight as it grazed the Sun. In both 
cases, GM/r(? <C 1, and the magnitude of the effects are indeed subtle. 

The advance of Mercury's perihelion is 574 arcseconds per century, only 43 of which are left 
unexplained by Newton's law after accounting for perturbations of Mercury's orbit caused by the 
other planets. In 1914, using the nearly completed GR theory, Einstein found that the rate of advance of 
the perihehon Aw « {2>G M / c^r orhit) '^orbit. which works out to 43 arcseconds per century. Einstein 
was thus able to account for a discrepancy that had been known for sixty years, and was dehghted, 
being moved to exclaim: 

For a few days, I was beside myself with joyous excitement. 

The quantitative confirmation of the deviation of a ray of Sun-grazing starlight brought worldwide 
acclaim to Einstein (an authoritative discussion can be found in [1]). In 1911, again before the final 
version of GR was in hand, Einstein predicted that the deflection should amount to 2GM j Rq(? = 0.87 
arcsec, where Rq is the radius of the Sun. This result could not have been especially satisfactory to 
Einstein, since it can be derived from Newtonian mechanics, if one assumes that light is corpuscular, 
composed of particles with a mass equivalent Ejc? . The result is referred to as the "Newton value," 
and is one half of the revised value derived by Einstein in 1915 (1.74 arcsec). In 1919, the result of 
a measurement of the deflection that confirmed Einstein's prediction was announced to the Prussian 
Academy, after which Einstein's reputation soared, both inside and outside academic circles, and the 
notions of warped space and time - warped spacetime - found their way into mainstream physics. 

The revolutionary implications of GR forced dramatic revisions of the concepts of space and time, 
not the least of which was the abandonment of the action-at-a-distance model of gravity. By contrast, 
the applications of GR in the early 1900s could be regarded as subtle corrections to Newtonian theory. 
However, discoveries of the 1960s — the discovery of compact X-ray binary sources [2], the recog- 
nition of the nature of quasars [3], and the discovery of radio pulsars [4] — provided examples of 
environments for which general relativity constituted substantially more than a gentle correction to 
Newtonian physics. 

The succession of rapid-fire astronomical discoveries of the 1960s marks the birth of relativistic 
astrophysics. Perhaps the most remarkable idea to present itself during that era was that black holes 
could exist in nature. In §2, a brief history of the black hole concept is presented. Here, we remark only 
that, in 1916, only a few months after Einstein published his final version of GR, Karl Schwarzschild 
discovered the solution to the field equations corresponding to the exterior of a sphericaUy symmetric 
mass distribution. This solution is the starting point for the study of black holes. Prior to the 1960s, 
however, references in the scientific literature to what eventuaUy came to be caUed black holes were 
sporadic and met with skepticism. 

Taken alone, GR does not prescribe a minimum black hole mass. However, those objects currently 
identified as black holes are believed to have formed from the catastrophic collapse of massive stars, 
and, combining GR with the equations of stellar structure, a minimum mass of about IMq is required 
[5] [6]. According to GR, a classical theory, matter undergoing gravitational collapse to form a black 
hole is crushed inexorably until it ends up in a point, known as the singularity. The singularity is 
enshrouded by a mathematical surface known as the event horizon, within which matter and light are 
trapped, destined to merge with the singularity. The event horizon has a radius 2GM/ (? (about 3 km 
per solar mass), if the black hole is not spinning. The event horizon can be as small as GMjc?, in the 
case of a black hole spinning at its maximum rate. A black hole is a simple object, in the sense that 
spacetime outside its event horizon can be described entirely by its mass, spin, and charge. Therefore, 
black holes belong to a three-parameter family - absurdly simple compared to, say, main-sequence 
stars. And yet, black holes so violate one's physical intuition that even Einstein doubted that they could 
exist in nature. In the words of Igor Novikov [7] : 
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Of all conceptions of the human mind perhaps the most fantastic is the black hole. Black 
holes are neither bodies nor radiation. They are clots of gravity. 

Outside the event horizon, the behavior of matter is dictated, in part, by the spacetime geometry. 
If radiating matter of sufficient quantity exists near the event horizon of a black hole, such that it 
is observable at Earth, then some of the effects of warped spacetime are observable, and we gain 
experimental access to physical environments for which GMjrc? is not negligibly small. This regime 
is known as the strong field limit, where relativistic effects are of fundamental importance. Nature has 
obhged us, by providing (at least) two classes of black hole systems - black hole X-ray binaries and 
active galactic nuclei (AGN) - in which prodigious quantities of matter are flowing "onto" the black 
hole, resulting in the release of copious electromagnetic radiation. In fact, detections of X-ray line 
radiation from regions where GMj r(? « 1 have been claimed. 

The flow of material referred to above is called accretion, defined as the capture of matter by an 
object's gravitational field, where it is presumed that the fate of the captured material is coalescence 
with the gravitating body, or, in the black hole case, passage through the event horizon. The physics 
of accretion has been an active area of research for over three decades. For the purposes of this paper, 
the importance of accretion is tied primarily to the role it plays in the conversion of gravitational po- 
tential energy into radiation. Owing to the prevalence of angular momentum in the cosmos, accretion 
often involves a disk. Accretion disks provide the means by which to dissipate angular momentum, 
allowing accretion, which is accompanied by the release of energy. The extraction of energy is espe- 
cially efficient if the inner edge of the accretion disk extends to small radii, so that the ratio Mj iiinner 
is large. This is the case for accretion onto compact objects - white dwarfs and, especially, neutron 
stars and black holes. The luminosities of the accreting compact objects that populate the Milky Way 
Galaxy range as high as ~ 10^ L©', while AGN can exceed that by seven, eight, even nine orders of 
magnitude [9]. 

Accretion disks present a number of theoretical challenges. The complexity of the problem, if one 
hopes eventually to work from first principles, is suggested by S. Shapiro and S. Teukolsky [10]: 

In the general accretion case, one must solve the time-dependent, multi-dimensional, rel- 
ativistic, magneto-hydrodynamic equations with coupled radiative transfer!" 

Given the context of the earher parts of this section, the above quote begs the question as to whether or 

not the strong field regime of GR can be tested through studies of accreting black holes. Considering 
that accretion onto black holes is one of the dominant light-producing mechanisms in the Universe, 
developing reliable working models of disk accretion is, from the astrophysics standpoint, hkely to be 
considered a higher priority than testing GR, although the latter is certainly a long-term goal. Moreover, 
it is probably fair to say that, at present, discerning possible inadequacies in strong field GR, a "clean" 
theory, using "messy" accretion disk models is problematic. In the meantime, recognizing that GR has 
stood up to every test to which it has been subjected, the theory in full is generally taken as being 
correct, and incorporated into the array of theoretical tools used to study black hole accretion. 

Accretion-powered objects are bright X-ray sources [2] [11]. The observed variability of the X- 
ray flux [12] [13] [14], rapid compared to variability at longer wavelengths, implies that X rays are 
produced preferentially in the inner regions of the disk. ^ The X-ray continuum spectra vary from 
source to source, and can vary in time as well, but for our purposes it is adequate to think of accreting 
black holes as exhibiting power-law continua that extend up to hard X-ray energies (~ hundreds of 



The symbol refers to the standard solar value of the subscripted quantity. Throughout the paper, we refer to the mass, 
radius, and bolometric luminosity of the Sun - Mq = 1.99 x 10^^ g; Rq = 6.96 x 10^° cm; Lq = 3.9 x 10^^ erg s~i 
[8]. 

^ The time scale for appreciable variations in the observed flux is associated with the Ught-crossing time of the source, so that 
the linear extent of the soiffce is given by R ~ ctcross- 
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Fig. 1. Broad iron Ka emission from two Seyfert 1 galaxies: MCG-6-30-15 {left panel) [26] and NGC 3516 
(right panel) [126] as recorded with the ASCA Solid-State Imaging Spectrometer (crosses) following subtraction 
of a continuum model. The line flux is scaled in multiples of 10~* photons cm~^ keV~^. Line profiles consist 
of a narrow emission component at 6.4 keV, the approximate rest energy of near-neutral iron ions, and a broad 
redshifted wing, believed to result from relativistic effects within a few gravitational radii of a supermassive black 
hole. Figures courtesy of K. Nandra (see [20]). 



keV) [15] [16]. If an optically thick accretion disk surrounds the compact object, and if the X-ray 
continuum flux irradiates the disk, then spectral signatures of the interaction of the X-ray continuum 
and the disk are expected to imprint themselves on the overall spectrum [17] [18] [19]. This imprinting 
is known as "reflection." Among the expected reflection features are a blend of X-ray lines from iron 
lying in the 6.4-7.0 keV spectral range. These lines are produced as the result of Is photoionization 
of multi-electron iron ions. It was also shown that, if the iron lines are generated very near a black 
hole, relativistic effects will skew the line, i.e., relativistic Doppler effects would broaden the line to 
several x 10,000 km s^^, and the gravitational redshift would produce observed emission down to 
several hundred to a few thousand eV below the rest energy (see Fig. 1), depending on the radial 
extent of the disk and its inclination with respect to the plane of the sky [21]. Precedence for the 
possibility of observing relativistic effects in emission lines was set by observations of the elliptical 
galaxy Arp 102B in the visible band, where the shapes of Ha and H/3 emission lines were analyzed in 
the context of a relativistic accretion disk model [22]. Broad X-ray iron line emission was detected in 
the Galactic black hole X-ray sources Cyg X-1 [23] and 4U1543-47 [24], possibly implying relativistic 
broadening, although the poor spectral resolution available at that time made definitive conclusions 
problematic. Alternative mechanisms for producing broad lines were shown to be unlikely [25]. The 
breakthrough observation was that of the Seyfert 1 galaxy MCG-6-30-15, using the ASCA observatory, 
which exhibited evidence for highly-broadened and skewed iron line emission [26] [27]. 

The discovery of skewed iron line emission has motivated a great deal of theoretical research, and 
has intensified observational searches for more and better X-ray data of black hole systems. By the 
current paradigm of black hole accretion disks [28], the study of iron line emission from the inner 
orbits of the accretion disk surrounding black holes provides the most direct view of the region within 
a few gravitational radii of the central engine. Our intention is to provide a primer covering the three 
topics that unite this phenomenon: (1) black holes; (2) accretion theory and modeling; and (3) X-ray 
spectroscopy and spectral modeling. Having here presented a quick synopsis, in the following sections 
each is elaborated in turn. It is fair to note that the relativistic disk interpretation is not unanimously 
accepted. Not surprisingly, with far-reaching consequences at stake, alternative models for producing 
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the broadened line have been and are being investigated (e.g., [29], [30]). However, a comparison of 
the successes and failures of the various models would take us too far afield, and we focus here on the 
relativistic disk interpretation. 

We also omit from detailed discussion the subjects of aperiodic and quasi-periodic variabiUty in 
the X-ray emission from black hole systems, subjects that are extremely rich from a phenomenological 
standpoint, as evidenced by data collected with the Rossi X-ray Timing Explorer (RXTE). In black hole 
X-ray binary systems, the observed time scales for variability are of the same order as the dynamical 
time scales expected for accretion flows near the black hole. Therefore, the temporal behavior of the 
X-ray emission provides another powerful view of the relativistic dynamics associated with the inner 
accretion disk. Although studies of the timing properties of black hole systems have, for the most 
part, constituted a somewhat "orthogonal" approach relative to spectroscopic analysis, there is a great 
deal of potential in analyzing the X-ray behavior of accreting black holes in the spectral and temporal 
domains when measured simultaneously. For reviews, see [13] [31] [32]. 

We have endeavored to provide the reader with a sufficient number of references to begin an ex- 
ploration of the literature, but we make no claim of unbiased completeness. The collection of papers 
cited is a reflection of our own experience with the subject matter. We note that excellent reviews of 
this topic have recently appeared [33] [34]. We take a more tutorial approach in this paper, and hope 
that readers will find our paper to be complementary, although there is some unavoidable overlap. 

The organization of the paper is as follows: In §2 we briefly sketch the historical development of the 
black hole concept, from Michell's dark star to Schwarzschild's solution of Einstein's field equations to 
the golden age of theoretical black hole research. Properties of the two size classes of black holes that 
are convincingly estabUshed - stellar-sized and supermassive - are outlined in §3. Also in that section 
we comment on the more speculative intermediate-mass black holes, and the microquasar phenomenon. 
A few of the core concepts associated with accretion power are treated in §4, including the a-disk model 
and the concept of accretion efficiency. In our treatment of general relativity, we forego a treatment 
of the theory proper, and instead focus on a few of the applications relevant to line spectroscopy of 
accreting black holes, including the implications of both the Schwarzschild (§5) and Kerr (§6) solutions. 
These include the gravitational redshift, the nature of circular orbits, and Ught motion. Vertical structure 
calculations of accretion disk atmospheres and the basic atomic physics processes that give rise to 
iron K spectra are the subjects of §7. Topics here include X-ray photoionization codes, the X-ray 
fluorescence mechanism, the fluorescence yield, resonant Auger destruction, and X-ray reflection. We 
conclude in §8. 

2. History of the Black Hole Concept 

Given the association of black holes with Einstein's GR, it might be natural to conclude that the 
notion of objects from which light cannot escape was engendered by 20th Century physics. However, 
allusions to the idea were already around hundreds of years earlier. In 1676, Olaus Roemer, noting 
variations in the period of Jupiter's moon lo, discovered that the speed of Ught is finite. In 1728, James 
Bradley, through observations of stellar aberration, produced further confirmation and a more accurate 
value of the speed of light, 295,000 km s^^. The concept of escape velocity, which, for a spherical 
mass M of radius R, is v^sc = {2GM/Ry^^, was present in Newton's time. In 1783, the English 
geologist John Michell combined these previous pieces of information, and realized that it would be 
theoretically possible for gravity to be so overwhelmingly strong that nothing - not even light - could 
escape. Michell's statements were made in the context of a corpuscular theory of light, the in-vogue 
approach at the time: 

Suppose the particles of hght to be attracted in the same maimer as all other bodies... [then] 
there should exist in nature bodies from which light could not arrive at us. 
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Michell went further, and proposed that from the motion of companion stars^ 

... we might still perhaps infer the existence of the central objects with some degree of 
probabiUty. 

Michell's objects were dubbed dark stars. Pierre Simon Laplace pubUshed similar comments in his 
1795 Exposition du Systeme du Monde, and added that 

... it is therefore possible that the greatest luminous bodies in the Universe are on this 
account invisible. 

Late 19th Century experiments fueled the rise of quantum mechanics and the wave interpretation of 
light, and the ideas of dark stars were soon forgotten. 

The mathematical discovery that would lead to the modern concept of the black hole was introduced 
in 1916, as the first solution to Einstein's field equations. At the outbreak of World War I, in August 
1914, Karl Schwarzschild volunteered for military service. He served in Belgium as commander of 
a weather station, in an artillery unit in France calculating missile trajectories, and in Russia. During 
this time he wrote two papers on Einstein's relativity theory (and one on Planck's quantum theory, 
explaining the splitting of the spectral lines of hydrogen by an electric field - the Stark effect). Einstein 
was later to say of Schwarzscluld's work: 

I had not expected that one could formulate the exact solution of the problem in such a 
simple way. 

Regrettably, while in Russia, Schwarzschild contracted a fatal disease and went home to die at the 
age of 42, never believing in the physical nature of what his solution implied for dark stars. In fact, a 
back-of-the-envelope calculation shows that for an object to exist within a sphere of radius 2GM/c^, 
its mean density must satisfy p > 3c®/327rG^M^, or p > 10^^ g cm~^ for a stellar-mass object, which 
must have seemed absurd at that time. 

With the end of the First World War, further astronomical tests of general relativity led to popular 
and scientific interest in Einstein's ideas. Most notable was an expedition led by Arthur Eddington, 
the goal of which was to obtain a measurement of the bending of starlight by the Sun during a solar 
eclipse. Meanwhile, research on stellar structure and evolution thrived during these years. In 1930, S. 
Chandrasekhar computed the first white dwarf models, taking into account special relativistic effects in 
the degenerate electron equation of state, based on the Fermi-Dirac statistics introduced in 1926, and 
on the ideas of W. Fowler regarding the equilibrium between electron degeneracy pressure and gravity. 
Chandrasekhar discovered that no white dwarf could sustain a mass larger than 1 AMq . It was apparent 
that if the electron degeneracy cannot withstand gravity, further collapse is inevitable. 

In 1934, W. Baade and F. Zwicky predicted that this further collapse would strip the atoms of their 
electrons, packing the nuclei together, while forming a neutron star. These stars were expected to be 
about 10 km in diameter, but with densities on the order of a billion tons per cubic inch. The neutron 
had been discovered by Chadwick only two years earlier. In this same work, Baade and Zwicky coined 
the term supernova and introduced its concept. In acknowledgement of one of the most prescient papers 
ever written, even the comic strip of the Los Angeles Times would conunent on the issue. Zwicky would 
later say of this insert: 

This, in all modesty, 1 claim to be one of the most concise predictions ever made in science. 
More than thirty years were to pass before the statement was proved to be true in every 
respect. 



' Michell is credited with the idea of binary stars. 
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In 1939, Oppenheimer and VoLkoff [35], and, a few months later, Oppenheimer and Snyder [36] 
realized that this inevitable stellar collapse implies a stellar evolutionary scenario for the formation 
of black holes. In the first paper, the first detailed calculations of the structure of neutron stars were 
performed, establishing the foundation of the general relativistic theory of stellar structure. The second 
paper focused on the collapse of a homogeneous sphere of a pressure-free fluid using general relativ- 
ity. Oppenheimer and Snyder's star was precisely spherical, non-spinning, non-radiating, uniform in 
density, and with no internal pressure. It was this set of assumptions, perhaps, that gave rise to skepti- 
cism regarding the applicability of the results to real imploding stars. The most impressive prediction 
in the Oppenheimer and Snyder work was the fact that an external observer - one far from the star 

- would see the implosion come to a halt at the event horizon, whereas one riding with the falling 
material would witness the whole process to the end, finding herself drawn into a point of infinite 
density, notwithstanding the meager hopes of surviving the trip until contact with the singularity. They 
concluded: 

When all thermonuclear sources of energy are exhausted, a sufficiently heavy star will 
collapse. Unless [something can somehow] reduce the star's mass to the order of that of 
the sun, this contraction wiU continue indefinitely. 

These discoveries in the field of stellar evolution directed attention to the Schwarzschild model 
of the spacetime exterior to a star. But yet again, a new world war would soon divert the efforts of 
most scientists towards military and nuclear research. Work on black hole physics recommenced in 
the mid-1950s, when D. Finkelstein discovered an alternative coordinate system for the Schwarzschild 
geometry that helped to clarify the Oppenheimer and Snyder results as being caused by time dilation 
in a gravitational field. J. A. Wheeler, first an opponent of the black hole concept, began to accept it, 
and tried to construct a quantum mechanical view of the inner singularities. FinaUy, Kerr discovered 
the spiiming black hole solution to the Einstein field equations in 1963 [148]. 

The term "black hole" was introduced in 1967 by Wheeler. In a 2003 interview with P. Davies, he 
explains his choice: 

The occasion was a meeting in the fall of 1967 at the Institute of Space Studies in New 
York to consider this marvelous work of Jocelyn Bell and Anthony Hewish on the pul- 
sars. What could be the cause of these absolutely regular pulses from some object out in 
space, and one obvious possibility was vibration of a white dwarf star, another was the 
rotation of a neutron star. But then 1 thought that to keep one's mind open, to look at all 
the possibilities one ought really to look at the gravitationally completely collapsed object. 
Well, the very words sounded so foggy, so ethereal, so far from touchability that nobody 
resonated to that as a possibility to be investigated. So in desperation I adopted the words 
Black Hole. Well, here at least was a name. 

The golden age of theoretical research in black hole physics started in 1964, and was to last for 
at least a decade. In this decade, computer codes used in hydrogen bomb research were adapted to 
study stellar collapse, topological methods and thermodynamical ideas were introduced into the study 
of black holes, the "no hair theorem" was demonstrated (a black hole has no characteristics indicative 
of the star from which it came, and that only three parameters - mass, angular momentum, and charge 

- are needed to describe it), and the cosmic censorship conjecture was stated (there are no naked sin- 
gularities, but rather, they are dressed with horizons), among many other results. After 1975, whereas 
important theoretical progress was still being made, the field slowly started to be dominated by an 
astronomical search for black holes in the Universe, at all scales. By that time, there were ample astro- 
physical reasons to expect objects described by solutions to Einstein's field equations to become part 
of the ontology of the Universe. 
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3. Accretion Power 

Accretion is defined as the capture of matter by an object's gravitational field, where it is presumed 
that the fate of the captured material is coalescence with the gravitating body, the results of which are 
an increase in the object's mass, as well as the conversion of gravitational potential energy into other 
forms. The introduction of accretion flows into astrophysics is often traced to the 1939 paper by Hoyle 
and Lyttleton [37]. Interestingly, its title is "The Effect of Interstellar Matter on CUmatic Variations," 
which examined the possibility that variable accretion of interstellar clouds by the Sun could be linked 
to the Ice Epochs of Earth's deep geological past. The notion of accreting black holes lay more than 
two decades in the future. The astronomical discoveries of the 1960s provided impetus for the rapid 
development of accretion theory, currently an active area of astrophysical research. We provide here 
a brief sketch of the theory as it pertains to our subject matter. A collection of seminal papers on the 
topic can be found in [38]. For a comprehensive textbook introduction to the physics of accretion, see 
[39]. An up-to-date review article, with a pedagogical approach, can be found in [40]. 



3.1. Basic Concepts 

To impart the fundamentals of accretion physics, we start with the Newtonian view of gravitation 
and a neutron star onto which matter falls. A test mass m that free falls from rest at "infinity" and 
comes to rest on the surface of a star of mass M and radius R, loses an energy GMm/R, which, for 
the purposes of this discussion, we assume to be converted into radiation. Therefore, given a more- 
or-less steady infall of matter, the energy generation rate — the accretion luminosity — is Lace = 
GMM^ccI R, where Mace is the mass accretion rate, the dot denoting differentiation with respect to 
time. It is sometimes convenient to express Lace in a special relativistic context, treating the energy 
release per unit time as being equivalent to a conversion of rest mass energy into radiation per unit 
time. Thus Lace = V^'^accC^, where rj is the dimensionless accretion efficiency. For the case of a 
spherical star, as discussed here, 77 = GM/Rc^. The appearance of relativistic correctness implied by 
the adornment of these last two expressions with factors of is entirely artificial, since the derivation 
of the accretion luminosity was obtained with purely Newtonian physics. Later, however, we will see 
that expressing Lace in terms of 77 follows quite naturally in the general relativistic context, and we 
continue to use 77 to relate Lace to Mace- 

Matter in an accretion flow is subject not just to gravitational forces, but to radiation forces, as well. 
The radiation produced near the surface of the accreting object exerts a pressure on the accretion flow. 
For radiation emerging from a gravitating source in a frequency interval dv, the force on a particle at 
distance r, in terms of an interaction cross-section (j{y), is 

where L^/^'Kr^hu is the areal rate at which photons in the frequency range [v, v -\- dv] arrive at r, and 
hv/c is the photon momentum. Thomson scattering sets the baseline for this integral, so that, ignoring 
relativistic corrections to the Thomson cross-section ctt, the net force, including gravity, is 

F = ^ - ^^"^P (2) 

assuming for simphcity a fully ionized hydrogen plasma (nip is the proton mass). For accretion to 
proceed, we must have F < 0, which imposes an upper limit to the accretion luminosity, known as the 

Eddington luminosity, which is given by 



(3) 
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By the same reasoning, there is an upper limit to the mass accretion rate. From the relation defining the 
proportionaUty between L and M, L = rjMc^, we define the Eddington accretion rate, 

Me = ^ = 1.4 X IQis M- — g s-\ (4) 

TjcaT -™0 

If a black hole accretes matter at the Eddington rate, then its growth rate dMsH/dt oc Mbh, and the 
mass grows exponentially with time, with an e-folding time Tgiowth — rjccTT / ^■nGnip ^ 5 x 10 yr, if 
r] = 0.1. 

Perhaps the simplest model of mass transfer in an X-ray binary system, to the extent that it is 
straightforward to obtain order-of-magnitude estimates, is that pertaining to a high-mass X-ray binary 
(HMXB). A more detailed description can be found in [41]. In HMXBs, a high-mass (O or B) star, 
radiates a bright UV continuum that transfers momentum to the stellar atmosphere through absorption 
in resonance lines, thus driving a mass outflow, or a stellar wind. In close proximity to the OB star is a 
neutron star, which captures part of the wind. In order to obtain numerical values of various quantities, 
we assume here that the mass of the neutron star is 1.4 Mq, and that its radius is 10^ cm, so that 
r] = 0.21. A simple estimate of the rate at which mass is captured by the neutron star can be found 
by finding the distance race at which the wind kinetic energy per unit mass, determined by the wind 
velocity Uw, is equal to the gravitational potential associated with the neutron star. This gives the 
accretion radius, r^cc = ^GM/ w^. Since the neutron star accretes at the rate Mace = 7rr^(.(.pw^'w, we 
can determine the accretion luminosity once we know the wind mass density at the position of the 
neutron star. This can be estimated by a mass continuity equation, assuming that the wind is spherically 
synmietric with respect to the OB star: = Atto^p^v^, where Mw is the total mass loss rate of the 
star, and a is the separation between the centers of mass of the components of the binary. Putting aU 
this together, we find 

f CM\ 

Mace ~ — Mw = 2 X 10^6 (^^)-4 g g-1 

To obtain the numerical estimate, we have written the various quantities in terms of multiples of typical 
values for those quantitities: a\2 is the binary separation expressed as a/ (10^^ cm); {v^)% is the wind 
velocity expressed as w/(10* cm s""'^); (Mw)-6 is the mass loss rate expressed as M^/ (lO^^M© yr~^). 

In obtaining an estimate of the mass accretion rate in a typical HMXB, there is no reason we could 
not have substituted a black hole for the neutron star. But let us now proceed to determine the resulting 
luminosity and radiation temperature. Using 77 = 0.21, derived above for the neutron star case, the rate 
of mass accretion found in Eq. (5) results in a luminosity Lace 10'^^ erg s~^, or approximately 10^ 
Lq. An order-of-magnitude estimate (actually a lower bound) of the radiation temperature associated 
with the accretion luminosity can be obtained by assuming that the radiant energy takes the form of a 
blackbody. If the energy is released uniformly over a spherical surface of radius R, which we assume 
to coincide with the neutron star surface, then Lace = 4:'jTR'^aT^^^, where a is the Stefan-Boltzmann 
constant. Thus we findTbb = {GMMg^^/^iraR^Y/^, or, again in terms of typical system parameters, 

n^ = i.2xl(f {v^)-^a-^'\M^fll K, (6) 

which implies that the emission peaks in the X-ray band. 

The simple model described above is consistent with observations of typical HMXBs. Is this model, 
on the other hand, relevant to black hole systems? For example, is it valid to simply evaluate 77 for a 
black hole by inserting, say, R = IGMjf? into t] = GM/Rc^, i.e., t/bh = 0.5? Clearly, the answer 
must be "no," since the calculation of 77 for a neutron star is based upon the assumption that the matter 
comes to rest on a material surface. In the cold fluid approximation, as used above, we might expect all 
of the energy release in black hole accretion to occur below the event horizon, leaving 77 = 0. In a more 
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realistic treatment, one assumes a finite gas temperature, and includes thermal energy balance in the 
calculations. Some fraction of the accretion energy is radiated away before crossing the event horizon. 
Such a scenario appUes to accretion of the interstellar medium by an isolated black hole. Detailed 
calculations [42] show, however, that for an ambient ISM density and temperature of 1 cm^'' and 10^ 
K, respectively, the efficiency rj ~ 10~^°, and the expected luminosity is expected to be only 10^^ 
erg s~^ per solar mass of the accretor, far too dim to be observed (cf., [43] [44]). The inefficiency of 
spherical accretion by a black hole is borne out by the realization that, in spite of the large number of 
black holes postulated to populate the Milky Way Galaxy, only the merest fraction (perhaps 1 in 10'') 
has been identified.'' 

In acknowledging that radial accretion onto a black hole is unUkely to produce a light source suf- 
ficiently bright that it is observable at Earth, we turn instead to accretion via a disk. We will see that 
disk accretion can be an extraordinarily efficient means by which to extract energy from matter near a 
black hole. 



3.2. Accretion Disks 

Angular momentum favors the formation of disks around accreting masses. The surfaces of con- 
stant gravitational potential between two stars in a binary system are such that a saddle point forms 
on the line of centers (the imaginary line connecting the centers of mass) between the two stars. The 
equipotential surface that intersects this saddle point is called the Roche lobe. If matter from one star 
reaches this saddle point, known as the inner Lagrangian point, or L 1 , then it can accrete onto the other 
star. Since there is relative motion between the stars, however, the overflowing matter cannot strike the 
accreting star directly, but rather goes into an orbit [45] dictated by its specific angular momentum. If 
sufficiently large quantities of matter spill through LI, however, these orbits cannot persist. Instead, 
collisions result in loss of energy and circularization of the orbits. Viscosity spreads the disk, transport- 
ing angular momentum outward and matter inward. This model describes the mass transfer mechanism 
in X-ray binaries, a pair of stars in close orbit, one of which is a gravitationally collapsed object — a 
neutron star or a black hole. The mass donor star is known as the companion. The X-ray binary phase 
of a Roche lobe overflow system is initiated when the companion expands to fill its Roche lobe, a con- 
sequence of normal stellar evolution [46]. The accretion of matter, if the overflow rate is sufficiently 
large, can power a highly luminous X-ray source. In this case, however, as contrasted to the HMXB 
case, the disk itself radiates much of the luminosity, irrespective of the presence of a material surface 
at the gravitating center. 

Disk accretion by Roche lobe overflow was first invoked to explain cataclysmic variable systems 
[47], in which a normal star transfers mass to a white dwarf (i?wd ^ 10^ cm) in a close binary system. 
The fundamental theoretical aspects of accretion disks in X-ray binaries were worked out by Pringle 
and Rees [48] and by Shakura and Sunyaev [49]. The latter authors derived an analytic model of disk 
flow, where the viscously dissipated energy is assumed to be locaUy radiated as blackbody emission, 
with temperature given by 



Tbb(r) 



3GMM 



1/4 



(7) 



where i?* is the radius of the accretor, or, in the black hole case, is the innermost stable circular 
orbit of disk material (see §6). Equation (7) is similar to the expression derived for HMXBs, showing 
that, for similar parameters. X-ray-emitting temperatures are expected. The energy dissipated per unit 

By "identified," we mean "inferred." The designation black hole candidate is applied when alternative explanations do not 
suffice. 
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area is 

from which Eq. (7) follows. The disk luminosity idisk can be found by integrating Q{r) over both 
disk faces, from i?» to some outer radius i?outer- A good approximation for Ldisk is found by letting 
-Router oo, yielding the result -Ldisk = GMM /2Rit [39]. If the accretor has a solid surface, the 
remainder of the energy is dissipated as it comes to rest on the surface. If the accretor is a black hole, 
then the remainder simply disappears behind the event horizon. In this latter case, if the black hole 
does not spin, the inner disk radius is taken to be &GM/c^ (see §5), which gives -Lace = M(? jVl. i.e., 
T] » 0.083. For a black hole spinning at its theoretical maximum rate, we set -R* = GM/c^, which 
gives 7? = 0.5. In §5 and §6, we derive relativistically correct values for the accretion efficiencies of 
non-spiiming and spinning black holes, respectively. 

In terms of accretion disk dynamics, the essential breakthrough introduced by Shakura and Sun- 
yaev [49] was the a-prescription, a parameterization of the vertically-averaged stress that assumes it 
can be written as aP, where P is the vertically-averaged total pressure in the gas (gas pressure plus 
radiation pressure), and where a is a constant. The a-prescription permits the solution of the coupled 
set of eight equations that determine the disk structure. Although the a-prescription sidesteps a good 
bit of complex physics, it turns out that many results are fairly insensitive to the value of a. For exam- 
ple, the disk central temperature, i.e., T(r, z = 0), is proportional to a~^/^, whereas, Eqs. (7) and (8) 
show that some results of interest have no dependence at all on a. For this reason, the a-disk has seen 
wide use. Indeed, the optically thick a-disk model has been accepted with only a few modifications 
since its inception in 1973 ([50] [51]; [52]). Among these modifications are the development of vari- 
ous branches of self-consistent accretion flow solutions such as advection-dominated accretion flows 
[53] [54] that successfully describe the spectral behavior of black hole XRBs [55] [56], advection- 
dominated inflow/outflow systems [57], and convection-dominated accretion flows [58]. Comparisons 
of the a-disk model with observations of X-ray binaries and AGN are presented in [39] and [59] [60], 
respectively. Modem accretion flow models are discussed in [61]. 

In spite of the successes and ease of use of the a-disk model, a physical understanding of a (or 
something like it) has been an important goal for accretion disk modelers, since it determines the es- 
sential mechanism that drives the transport of angular momentum. The viscous mechanism in accretion 
disks allows accretion to occur by transporting angular momentum outward, matter inwards, and by 
dissipating gravitational energy into heat inside the disk. The viscosity mechanism has been identified 
in theory as a magneto-rotational instability arising from the entanglement of magnetic fields caused 
by differential rotation of gas in Keplerian orbits [62]. In MHD models, this mechanism provides the 
energy dissipation and angular momentum transfer needed to naturally produce mass accretion with 
self-sustained magnetic fields B that are much smaller than the equipartition level (B'^/Stt <C pv"^), 
where p is the gas mass density and v is the thermal velocity. Numerical magnetohydrodynamic (MHD) 
models [63] show that a is not constant, but that it ranges from 10^^-10"^. 

While the X-ray continuum in accreting neutron star systems is interpreted, in part, as thermal 
emission from the disk, the presence of (apparently) non-thermal continuum radiation in these sources, 
and AGN, as well, is not naturally explained with the a-disk formalism. In fact, the high-energy emis- 
sion is especially problematic in AGN, if we consider Eq. (7). Anticipating a result that is derived in 
§6, set -R* = QGM/c^, which corrresponds to the irmer disk edge in a non-spitming black hole, and 
express the mass accretion rate as a fraction /e of the Eddinglon accretion rate, M = JeMe- Finally, 
if we express r as a multiple of the gravitational radius GM/ (?, so that x = (?r/GM, then Eq. (7) can 
be written as 



7bb(a;) 



2r]aaTGM 



(9) 
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The key result to take away from this expression is the scaling Tbb(a;) oc M~^/**. Thus if we accept 
the a-disk model for stellar-sized accreting black holes (Mbh ~ WMq), we are faced with a glaring 
inconsistency if we wish also to apply it to AGN (Mbh ~ W^Mq); all other things being equal, 
^br/^bb' = (M^h/M^b)^^*' which means that a 10^ K blackbody for an X-ray binary is scaled 
down to about 3 x 10^ K for an AGN disk. The AGN emission is thus expected to peak in the UV 
band, with very little luminosity appearing in the hard X-ray region. While the UV peak, the "big 
blue bump," is a well known component of AGN spectra [64] [65], the a-disk model says nothing 
about the X-ray flux that characterizes AGN [66], and refinement is required in order to match disk 
models to observations. For example, a more realistic accounting of radiation transport effects in the 
disk atmosphere partially ameliorates this problem [67] [68] [69]. 

The favored explanation of the X-ray luminosity of accreting black holes posits that the blackbody 
component originates in a cold, optically thick accretion disk, whereas the hard X-ray power-law com- 
ponent is produced in an optically thin hot corona by thermal Comptonization of disk photons [70] [71]. 
Comptonization refers to the deformation of a radiation field as it interacts through Compton scattering 
with an electron distribution, where a self-consistent solution determines both the spectral shape and 
the electron temperature [72] [73] [74]. Constructing theoretical models of continuum production is 
complicated by the need to include 6+- e~ pair production in a manner that is self-consistent with the 
radiation field [75] [76] [77] [78]. 

In disk accreting systems where a hard X-ray source is present, the disk is exposed to this radiation 
and will be heated by it. In fact, radiative heating can exceed internal viscous heating in some regions 
of the disk. The temperature structure of the disk can thus be controlled by the X-ray field, photoioniz- 
ing the gas, suppressing convection, and increasing the scale height of the disk [79]. Photoionization in 
the disk is balanced by radiative and dielectronic recombination [80], and possibly three-body recom- 
bination [81] and charge transfer recombination [82]. These recombination processes produce discrete 
Une emission, which constitute spectral components distinct from the fluorescence component of the 
spectrum, potentially providing corollary information relating to the disk structure. A quantity that is 
commonly used to describe emission lines that are superimposed on a continuum is the equivalent 
width. Since we make several references to the equivalent width in later sections, we define it here for 
convenience. The experimentally measured line equivalent width is defined by 



The equivalent width is usually quoted in units of eV or keV. Theoretical predictions of follow 
from calculations of accretion disk models, which can then be compared to observations. For example, 
early model calculations of the expected iron Ka equivalent width, based on the response of an X-ray 
irradiated slab, showed that 90-150 eV was attainable [83]. 

The production of line emission, regardless of the formation mechanism, and the transfer of these 
lines through the overlying atmosphere requires, in principle, a detailed calculation of the opacity 
distribution along their lines of flight. The spatial distribution of the line emissivity, as well as a proper 
accounting of the probability of escape must, therefore, account for the vertical structure of the disk. 
This includes, in the Umit that the vertical structure is in steady-state, the density, temperature, and 
charge state distribution. Thus each annulus of the disk can be thought of as a stellar atmosphere, 
where the equations of energy flow must be solved self-consistentiy with the ionization equations, 
possibly constrained by hydrostatic equilibrium. The explicit vertical structure of the disk, however, 
does not follow directly from the a-disk model. For example, the surface density - the line integral of 
the density in the vertical direction - appears as one of the variables in the disk equations. We return to 
the subject of the vertical structure of accretion disks as it bears on X-ray fluorescence in §7. 




(10) 
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4. Astrophysical Black Holes 

In this section, we present a brief survey of black holes from an astronomical point of view, touching 
on the gross characteristics of black holes in X-ray binaries and AGN. 

4.1. Black Hole X-ray Binaries 

Currently about 250 X-ray binaries are known in our galaxy [84] [85], possibly representing an 
underlying population of more than 1000 objects. X-ray binaries can be separated into two populations: 
the low-mass X-ray binaries (LMXBs), where "low-mass" refers to the companion star to the compact 
object, an older population, concentrated near the Galactic bulge; and the HMXBs (see §3.1), younger 
systems, concentrated in the spiral arms. Mass transfer from the companion to the compact object 
is found to typically differ for these two classes as well. In general, LMXBs transfer mass to their 
companions through Roche lobe overflow (Fig. 2). Mass transfer in HMXBs is generally mediated by 
stellar winds. 

A subset of X-ray binaries, twenty or so, are identified as black holes paired with a star that serves as 
a mass donor for accretion [55] [86]. The most famous example of this class is Cygnus X-1, identified 
in 1972 as a black hole candidate [87]. A more common subclass of X-ray binaries are those for which 
the compact object is a neutron star, rather than a black hole. When the compact object in an X-ray 
binary system is shown to be more massive than about 3Mq, the compact object is a good black hole 
candidate. For a small group of X-ray binaries, the mass measurement has been performed with high 
precision. This group of X-ray sources is usually referred to as dynamically confirmed black holes. ^ A 
continually updated Ust of such objects can be found, thanks to J. Orosz^ (see [32] for a more detailed 
discussion). Among the dynamically confirmed black hole X-ray binaries, there are several recurrent 
X-ray novae, all of them LMXBs, and only three persistent sources - Cyg X-1, LMC X-1, and LMC 
X-3 - the latter of which are HMXBs. These special binaries span a large range in the parameter space 
of basic properties; for example, XTE J1118H-480 has Poib = 0.17 days and a binary separation of 
w 2.8i?o, whereas GRS 1915H-105 has Porb = 33.5 days and a binary separation 95Rq. 

Clearly, if a less strict standard of evidence is accepted, many more black hole candidates can be 
identified (see, e.g., [32] for a recent discussion). Still, given the small number of black hole candidates, 
it may be surprising that, from stellar evolutionary considerations, together with the incidence of the X- 
ray binaries that are beUeved to contain black holes, it is thought that the Milky Way Galaxy contains 
about ~300 million stellar-mass black holes [88] [89] [90] [91]. This number implies, assuming ^ 
IOM0 per black hole [32], that about 4% of the total baryonic mass (i.e., stars plus gas) of the Galaxy 
is in the form of black holes. 

The X-ray spectra of X-ray binaries are dominated by continuum radiation. The X-ray luminosities 
of persistent XRBs are intrinsically variable, and range anywhere from (typically) ~ 10"^^ erg 
to as high as the Eddington luminosity of few x 10^* erg s~^. In black hole X-ray novae, the range 
in X-ray luminosities is more dynamic, with quiescent luminosities as low as ^ 10^" crgs~^ 
[92]. Based upon spectral observations in the soft X-ray band (say, 1-10 keV), two main categories 
of emission states have been distinguished, depending mostly on the slope of the power-law (Fe oc 
E~°') describing the continuum emission of the system. Here, the differential photon count rate is 
given in units of photons cm~^ s~^ keV~^. The so-called low state (relatively low photon count rate) 
features non-thermal X-ray flux, typically with slopes near 1.6 - 1.7. Because of the hardness (i.e.. 



' Therelevantexperimentally inferred quantity is the mass function, defined by /(M) = Poji,i<r|/27rG=Mi sin'^ 

where Porb is the orbital period and K2 is the semi-amplitude of the velocity curve of the secondary, Mi is the black hole 
mass, i is the orbital inclination angle, and q = M2/M1, where M2 is the mass of the secondary, is the mass ratio. The 
equation for f{M) impUes that the value of the mass function is the absolute minimum mass of the compact primary. A 
secure value of f{M) may be sufficient to show that the mass of the compact X-ray soitfce is at least 3Mq. 

* http://mintaka.sdsu.edu/faculty/orosz/web/ 



© NRC Canada 



14 



unknown Vol. 




Fig. 2. Current basic understanding of the physical components and sites of emission in an X-ray binary system. 
In this example, the companion loses mass to the compact object through Roche lobe overflow. The accreting 
matter forms a disk, which can be thought of as a series of annular rings in quasi-Keplerian orbits, whose radii 
decrease under the action of stresses between neighboring annuli. Radiant energy is released, giving rise to X-ray 
emission at small radii. Hard X-ray emission impinges on the optically thick disk, which "reflects" a fraction of 
the incident flux, with spectral imprinting in the form of lines, absorption, and deformed continuum. A fraction 
of the accretion energy is redirected into kinetic energy, giving rise to the jets. A hot corona sandwiches the disk 
at small radii, which can affect the reflection spectrum through Compton scattering. From Fender and Maccarone 
[130]. 



relatively high ratio of hard X-ray flux to soft X-ray flux) of such indices, this state is also referred 
to as the low/hard state. The high state, by contrast, is characterized by intense quasi-thermal flux. In 
this state, most of the radiated energy is concentrated in a blackbody component, while the power-law 
contribution becomes softer, with power-law indices typically larger than 2. The usual interpretation is 
that the blackbody component originates in a cold, optically thick accretion disk, whereas the power- 
law component is produced in an optically thin hot corona by thermal Comptonization of disk photons 
[71] [98] [99] [100]. 

The first broad Fe Ka line observed was reported in the spectrum of Cyg X-1, based upon EXOSAT 
data [23]. Recently, the Chandra X-ray observatory was used to observe Cyg X-1 with the High Energy 
Transmission Grating Spectrometer in an intermediate X-ray state [ 1 02] . A narrow Fe line was detected 
at = 6.415 ± 0.007 keV with an equivalent width (see Eq. 10) of Wt = 16^1 eV, along with a broad 
line at £' = 5.82 ± 0.07 keV with We = MO'^Iq eV. A smeared photoelectric edge was also detected 
at 7.3 ± 0.2 keV. These results are interpreted in terms of an accretion disk with irradiation of the inner 
disk producing the broad Fe Ka emission line and edge, and irradiation of the outer disk producing the 
narrow line. The broad line is thought to be shaped by Doppler and gravitational effects and, to a lesser 
extent, by Compton reflection (see §7). 

Early studies of relativistically smeared Fe Ka lines from X-ray binaries relied on proportional 
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counter detectors with relatively poor spectral resolution (e.g., the Ginga's LAC and the RXTE's PCA, 
have an energy resolution of » 1.2 keV at Fe Ka). The response matrices of the detectors are uncer- 
tain at the 1-2 % level, while the Fe line profile is typically only 1-5 % above the X-ray continuum 
[32]. Therefore, interpretations of results from these instruments must be approached with caution. 
BeppoSAX, with a resolution of « 0.6 keV at 6.4 keV, has also been used to observe several black hole 
candidates. The observed iron emission generally appear to be rather symmetric, and may be more a 
product of Compton scattering than relativistic broadening. However, in the cases of GRS 1915+105 
[103] and V4641 Sgr [104], the iron lines are skewed, possibly implying relativistic smearing. In addi- 
tion, an observation using the XMM EPIC-MOSl detector led to the report of a broad, skewed Fe Ka 
emission in XTE J1650-500, perhaps suggesting the presence of a rapidly spiiming black hole [105]. 

4.2. Active Galactic Nuclei 

Approximately one in a hundred galaxies shows evidence for energy output that appears to be 
unrelated to normal stellar processes - rapid bulk motions, a bright, non-stellar radiation continuum, 
rapid aperiodic variability, and high-luminosity emission lines originating near the galactic nucleus. 
Members of this subset of galaxies are known as active galaxies, and the nuclei are referred to as AGN 
[106]. The classification scheme for active galaxies is somewhat daunting. We list here the various 
classes as they have been designated in the literature [107]: (1) radio galaxies; (2) radio quasars; (3) BL 
Lac objects; (4) optically violent variables; (5) radio quiet quasars; (6) Seyfert 1 galaxies; (7) Seyfert 
2 galaxies; (8) low-ionization nuclear emission-line regions (LINERs); (9) nuclear H II regions; (10) 
starburst galaxies; and (11) strong IRAS galaxies. Not all among this group are active by virtue of 
black hole accretion. For example, starburst galaxies are sites of anomalously large star formation 
rates, probably the result of a gravitational encounter with another galaxy, leading to an enhanced IR 
luminosity. We do not delve here into the phenomenology that gives rise to this complex taxonomy. 
In what follows, we reserve the term AGN for active galaxies beheved to harbor a supermassive black 
hole. 

Within the unification scheme [108], the underlying model for all classes of AGN is intrinsically 
similar. At the very center of the galaxy sits a supermassive black hole (~ 10^ - 10^° Mq), which 
accretes galactic matter through an accretion disk. Broad emission lines are produced in clouds orbiting 
above the disc at high velocity (the Broad Line Region [109]), and this central region is surrounded 
by an extended, dusty, molecular torus. A hot electron corona sandwiches the inner regions of the 
disk, probably playing a dominant role in generating continuum X-ray emission. Two-sided jets of 
relativistic particles emanate perpendicular to the plane of the accretion disc, the generation of which 
is still not fully understood. 

As pointed out in [110], X-ray surveys are an excellent means by which to locate AGN. Surveys 
using the Chandra and XMM-Newton observatories find a remarkable ^ IQ-^ X-ray sources per square 
degree, the vast majority of which are undoubtedly AGN [111] [112]. To date, however, the most 
secure detections of supermassive black holes (accompanied by a "rule out" of alternative models, 
such as the existence of dense clusters of stars or exotic particles) come from stellar proper motion 
in the Galactic center and the H2O megamaser of the nearby Seyfert 2 galaxy NGC 4258 [113] (see 
[114] and references therein for further discussion). Optical stellar and gas dynamical studies, generally 
using the Hubble Space Telescope, have revealed a large concentration of nuclear mass in several tens 
of candidates, beheved to be black holes. It is true, however, that these methods currently lack sufficient 
angular resolution to probe the spacetime at distances on the scale of the black hole horizon (e.g., see 
[115]). 

It is worth noting that nuclear activity is not a prerequisite for the presence of a black hole. In fact, 
most of the galaxies probed for supermassive black holes are not really active, but dormant quasars, 
relatively close to Earth, and for which the Keplerian signatures could be more easily discovered be- 
cause of the higher spatial resolution [116]. The most famous example is provided by our own Galaxy 
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[117]. After charting the kinematics of stars swirling around the central regions of the Milky Way [118] 
[119] [120], it has been found that the total mass of the region enclosed within a radius of 2 x 10^^ 
cm is approximately 3.7 x 10^ Mq, far more compact than what is possible for a stable distribution of 
individual objects; total gravitational collapse is required, according to GR. Indeed, the only remaining 
alternative candidate, other than a black hole, to describe the behavior of the innermost stellar orbits in 
the center of the Milky Way Galaxy comes from particle physics, and is known as a boson star [121]. 
Mass estimates for central supermassive black holes in about twenty nearby galaxies are also available 
(see review in [122]). Therefore, if we assume that the Milky Way Galaxy is not unique in this respect, 
it is possible that a large fraction of galaxies, active or not, harbor supermassive black holes near their 
dynamical centers. 

Currently, the Seyfert 1 galaxies MCG-6-30-15 [123] [124], Mrk 766 [125], and NGC 3516 [126] 
provide the most robust detections of the relativistically broadened iron Ka line. In addition to the 
preceding papers, a concise set of case studies of these three objects is provided in [127]. 

4.3. ]VIicroquasars and Jets 

Some X-ray binaries are sources of jets, high-velocity streams of oppositely directed particles orig- 
inating near the compact object. These sources are called microquasars [128] [129], exploiting an 
analogy discussed in more detail below. Figure 2, from [130], is a sketch of an X-ray binary, present- 
ing the major physical components and sites of emission in such systems. The basic underlying idea 
behind the analogy between quasars and microquasars is that the physics in all black hole systems is 
essentially governed by scaling laws, whose order parameter is the black hole mass. For instance, the 
scales of length and time of black hole related phenomena are proportional to the mass of the black 
hole. For a given critical accretion rate, the bolometric luminosity and length of relativistic jets are also 
proportional to the mass of the black hole. For a black hole of mass M the density and mean temper- 
ature in the accretion flow scale with and Af"^/"*, respectively. The maximum magnetic field at 
a given radius in a radiation dominated accretion disk scales with M~^/^, which implies that in the 
vicinity of stellar-mass black holes the magnetic fields may be 10^ times stronger than that found near 
supermassive black holes [131]. 

One of the most impressive phenomena occurring in both quasars and microquasars is the ejection 
of blobs of plasma at apparently superluminal speed [132]. Very Long Baseline Interferometry at radio 
wavelengths allows position measurements down to milli-arcseconds, thereby permitting detections of 
small changes in position on the sky. The apparent velocity is obtained by multiplying the observed 
proper motion by the distance to the source. Such superluminal sources were discovered as radio- 
galaxies and quasars, where the central black hole is supposed to have millions of solar masses. It was 
found that this very same phenomenon could also occur in Galactic sources, after the report on GRS 
1915-H105 by Mirabel and Rodriguez [128]. 

4.4. Intermediate Mass Black Holes 

Although consensus has not yet been reached, there are two types of data suggesting the existence 
of intermediate mass black holes, with masses between ^ 2OM0 and several thousand Mq [136]. First, 
there are numerous X-ray point sources, dubbed ultmluminous X-ray sources, that are not associated 
with AGN, and that have fluxes far beyond the Eddington limit of a stellar size black hole system. 
Second, several globular clusters show clear evidence for an excess of dark mass in their cores, which 
appears to be a single object. The number of such intermediate objects, as well as their actual existence, 
is yet under debate, and strongly depends on the mechanism that forms them. For both intermediate 
and supermassive black holes, the formation processes are not as well understood as they are for their 
stellar counterparts. 
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5. Black Hole Accretion In the Schwarzschlld Metric 

In this section we present in some detail a few of the results of general relativity that pertain to 
black hole accretion. As mentioned earlier, there are only three vacuum solutions for black holes. 
The simplest - the Schwarzschlld metric - describes the geometry of spacetime outside the event 
horizon of a black hole with zero charge and angular momentum; mass alone describes the geometry. 
A clear and concise derivation of the Schwarzschlld solution can be found in, for example, [137]. 
Before proceeding to the Schwarzschild solution, we introduce the unit conventions that have been 
adopted for deaUng economically with calculations involving relativity. 



5.1. Geometrized Units 

In the domain of special relativity, space and time are linked by the invariant spacetime interval. 
By convention, we can rewrite the familiar expression for the proper time interval in fiat spacetime 
(ir^ = dt'^ — drx? as rfr^ = di^ — rfx^ by expressing time in cm; 1 cm of time is, in conventional 
units, the time required for light to travel 1 cm, or iconv = icm/(3 x 10^'' cm s~^). 

In general relativity, the dimension of length is given to other quantities as well. For example, 
as mentioned in §1, the characteristic length scale in black hole physics is the gravitational radius^ 
Rg = GMjc?. Thus one often encounters dimensionless terms of the form rj {GM/c^). In adopting 
geometrized units, we make the transformation 

r r 

(11) 



GM/c2 M' 

where M carries the unit cm, so that r/M is dimensionless. Therefore, to convert mass expressed in 
cm to mass expressed in grams, use 



c 



2 



M(g) = - M(,^) or M(g) = 1.347 x 10^« M(,^). (12) 

For example, the mass of the black hole near the center of the Milky Way Galaxy (M « 3.7 x 10^ Mq 
or 7.4 X 10^® g in conventional units) is 5.5 x 10^^ cm. 

Angular momentum, symbolized here as J or L, has dimensions of length squared in geometrized 
units. The conversion is 



G 



"^(cin2) = -^(gcm^s-i)- (13) 



For convenience with mathematical manipulations, the spin of a black hole is characterized by the 
quantity a = J/M, where J is the angular momentum of the hole. The conversion between geometric 
units and c.g.s. units is 

^ _ '^(cm^) ^ -^(gcm^s-i) ^^^^ 
^(cm) ^(g)C 

so that spin has dimensions of length. Since mass also has dimensions of length, it is often useful to 
define a dimensionless spin parameter, a* = a/M. If the mass and angular momentum of a spinning 
object are known in conventional units, a dimensionless a/M follows from 

__Cl___C_J_| 

M~ GM^ l-s- ^ ' 



^ As an easy-to-remember conversion, the gravitational radius of an object is 1.477 km per solar mass. 
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5.2. The Schwarzschild Metric 

A spacetime metric provides a "rule" by which to relate measurements of events in different frames. 

The Minkowski metric, 

dr^ = - {dx^f - {dx^f - (16) 

can be expressed in a more concise form if we define a set of metric coefficients 77^/3. Thus one can 
write dr^ = —rjajs dx^dx^, where dx° = dt, rjoo = —l,rijj = 1 for j = 1,2, 3, and where summation 
is performed over all repeated indices. 

More generally, a spacetime metric has the form 

dr"^ = -g^.^dx^'dx", (17) 

where the metric coefficients can be functions of the spacetime coordinates. The Schwarzschild solu- 
tion to the Einstein field equations yields the metric 

dr^ = (1-^) dt^ - (1-—) ' dr^ - {dO^ + sin^ 6 dd^% (18) 



r J \ r 

which is the vacuum solution to the spacetime outside of a spherical mass Al (for a derivation, see, 
e.g., [137J L138J). The metric coefficients g^u can be read off directly. The f-coordinate is the time as 
measured by a distant observer. The r-coordinate, called the reduced circumference, is defined such 
that the circumference of a circle centered on the gravitating mass is precisely 27rr. The 6- and (/)- 
coordinates are defined such that r d9 and r d(f>, respectively, measure differential distances along tan- 
gents to a circle at r, again, centered on the gravitating mass. For r ^ 2M, the Schwarzschild metric 
reduces to the Minkowski metric (Eq. 16), after a transformation between spherical polar coordinates 
and Cartesian coordinates. There is, however, no transformation of coordinates that can globally re- 
duce the Schwarzschild metric to the Minkowski metric. Thus the spacetime described by the former 
is said to be curved, whereas the Minkowski spacetime is said to be flat — "flat" in the sense that it is 
quasi-Euclidean. 

Since the metric is spherically symmetric, the motion of test particles and photons is restricted to 
a plane. This plane, with no loss of generality, can be chosen to be the equatorial plane (9 ~ 7r/2). 
For events occurring within the equatorial plane, we can also set d9 — in Eq. (18), which gives the 
simpUfied proper time interval, 



(19) 



..^=(l-^).,^-(l-^)"^.^-.^.^^ 
The proper length interval ds {ds^ = —dr'^) is given by 

ds={l- — \ dr. (20) 



5.3. Gravitational Time Dilation and Gravitational Red Shift 

Even before Einstein arrived at his final formulation of the field equations, he argued the case for 
gravitational time dilation. The standard version of his derivation appeared in 1911, although he was 
aware of the effect as early as 1907 [1]. We provide here a derivation that follows the one presented in 
[137], except that we will take advantage of the explicit form of the Schwarzschild solution, thereby 
deriving a result that will be useful later. 

Imagine that two identical clocks are manufactured in a region of flat spacetime and calibrated 
so that every hour the clocks chime together, each clock thus marking proper time intervals At 
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corresponding to one hour. Two experimenters, A and B, are enlisted. A moves one clock to a po- 
sition at distance from a gravitating object. B moves the other clock to a position at distance 
tb > ta from the mass (Fig. 3). At each chiming, A records the time, and will maintain that the 
passage of time from chime to chime At^ is precisely one hour. According to Eq. (19), A measures 
At^ = At/ — 2M/rA- Light pulses are sent from ta to rs at each chiming of ^'s clock, and B 
records the arrival times. Light pulses require time - according to observers in any frame - to travel 
from A to B. Nevertheless, B records a pulse-to-pulse time separation of precisely Af a, provided that 
the hght paths are identical for each pulse. When he compares the chime-to-chime record of his own 
clock, however, he finds that Ats < At^, where B will insist that Ats is one hour. B measures 
AtB = At/ -^1 — 2M/rB- B is forced to conclude that time runs slower at than at rs- If the 
experiment were altered such that B sent pulses to A, then A would be forced to conclude that time 
runs faster at rs than at r^. This effect is known as gravitational time dilation. 

For our purposes, it is more appropriate to express time dilation in terms of light frequencies. 
Referring to the experiment above, and replacing the clocks by radiating atoms, label by i/obs the 
frequency of hght originating at A as observed at B, and label by the frequency of an identical 
source at B as observed by B. The ratio of the frequencies is, therefore. 



^rec 



AtB _ ( l-2M/rA y'^ 

I-^M/vb) • ^ ^ 



I^em AtA 

The case of primary interest - observing radiation from a distant source - corresponds to Tb oo. 
This gives us the expression for the gravitational redshift for a Schwarzschild black hole, according to 
the observer at infinity, 

1 - — ) • (22) 

As r ^ 2M, the observer at infinity finds that z/obs — * 0. For this reason, the radius r = 2M is 
sometimes referred to as the surface of infinite redshift. It is important to remember that this expression 
is valid only for the case of a stationary source and a stationary receiver. Below, we will accommodate 
the case of a source in orbit around a black hole. 



5.4, Conserved Quantities 

Before discussing the motion of particles and photons in the Schwarzschild metric, we will need 
two conservation laws, which foUow from the equations of motion. The equations of motion of free 
particles in curved spacetime can be derived from the Euler-Lagrange equations, for an appropriate 

Lagrangian A, according to 

d dA dA _ 

dp di°' ~ 'dx" ~ ^ ^^^^ 

where the dot denotes differentiation with respect to a timelike parameter p. A useful form for the 
Lagrangian is 

1 dx" dx'^ 

with p = t/iti., as shown in [10]. For the Schwarzschild metric, again reducing the problem to motion 

in the 6* = 7r/2 plane, this is 
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time at point A 



time at point B 



>- 

increasing distance from gravity source 

Fig. 3. Gravitational time dilation. Two identically manufactured clocks are positioned at points A and B, with 
rs > ta with respect to a gravitating mass. An observer at A sends pulses to observer at B at each chiming of 
the clock at A. Observer at B notes that chime frequency is higher at B than at A and concludes that time runs 
slow at A. See text for futher discusion. 



Evaluating Eq. (23) for the t-coordinate, we find 

E ( 2M\ dt 

- = 1 — (26) 

TO \ r J dr 

where the constant of the motion, labeled E/m, is identified with the energy per unit mass of the 
particle, since Eq. (26) reduces to the special relativistic expression E — m dt/dr for large r, i.e., the 
time-component of the momentum four-vector. Similarly, evaluating Eq. (23) for the (^-coordinate, we 
find a second constant of the motion 

iL^r'^. (27) 
TO dr 

where L/m can be identified as the angular momentum per unit mass. 



5.5, The Effective Potential 

The usual Newtonian approach to the problem of calculating the orbit of a test particle under the 
influence of a spherically symmetric gravitational field involves finding two constants of the motion - 
energy and angular momentum - and recasting the radial equation in terms of an effective potential. 
The effective potential is the sum of the gravitational potential and a centrifugal term. From an analy- 
sis of the radial excursions in this effective potential, one finds the two basic orbit classes - elliptical 
and hyperbolic - corresponding to bound and unbound motion, respectively, with circular orbits cor- 
responding to a special case of elliptical orbits. A similar, and quite fruitful, approach is adopted to 
analyze the orbit classes for motion in the Schwarzschild metric. 

To derive the effective potential for the Schwarzschild metric, notice that Equations (19), (26), and 
(27) can be treated as a set of three equations in the four unknowns dr, dt, dr, and dcj). Therefore, we 
can eUminate dt and d(j) in Eq. (19) to obtain 

(dr) (to) r) 



{L/m) 



(28) 
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Fig. 4. (Left panel) Effective potential for a test particle trajectory in the Schwarzschild metric (Eq. 29) with 
L/mM = 3.95. Dashed lines correspond to five different energies, increasing from bottom to top. See text for 
a discussion. (Right panel) Effective potential for various values of the angular momentum: from bottom to top, 
L/mM = [0., 2.5, 2%/3, 3.7, 4.0, 4.3, 4.6, 4.9]. Intersections of the dashed line with potential curves marks 
positions of stable circular orbits. The innermost stable circular orbit occurs at r/M = 6, corresponding to 
L/mM = 2^3. 



and one defines the effective potential according to 



V\ _ f 2M 
ml \ r 



(29) 



Note that the effective potential cannot be associated with an actual potential energy, as in Newtonian 
mechanics; in general relativity, it is impossible to separate energy into a kinetic term and a poten- 
tial term. Rather, the effective potential merely allows us to characterize certain aspects of particle 
trajectories without integrating the equations of motion. 

One example of the effective potential is plotted in the left panel of Fig. 4; as r ^ oo, V{r) 1 
as can be seen from Eq. (29). The curve is, in fact, qualitatively similar to analogous curves derived 
from the Newtonian prescription, except for small radii, where the two differ dramatically, i.e., the 
repulsive centrifugal barrier that extends to arbitrarily small radii in the Newtonian potential is absent 
at small radii in the Schwarzschild effective potential. Any inward bound particle that finds itself within 
the region of this downturn of V, which has been referred to as "the pit in the potential" [170], is 
destined to fall through the event horizon. Energies corresponding to such capture orbits are indicated 
by the top two dashed lines in Fig. 4. The top line, with E/m > 1, is labeled "capture/escape orbit" 
because inward bound particles are captured, whereas outward bound particles can escape to infinity. 
The energy level labeled "capture orbit" coiTesponds to a particle that, if inward bound, passes through 
the event horizon, and, if outward bound, aiTives at apastron, returns, since E/m < 1, then passes 
through the event horizon. Below that is the level for a particle with energy corresponding to an unstable 
circular orbit at r = rug. Test particles with this energy, whether directed inward from r > r^s or 
directed outward from r < r^s, will end up on this orbit, and remain there until perturbed. The next 
orbit is bound, since particles with this energy transit between two turning points. The Newtonian 
analogy for this orbit is an ellipse. In the Schwarzschild metric, the azimuthal oscillation frequency 
exceeds the radial oscillation frequency, leading to an orbit with a precessing apastron, with prograde 
precession. The final trajectory class, corresponding to the local minimum of V, is a stable circular 
orbit, discussed in the next section. 
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5.6. Circular Orbits in the Schwarzschild Metric 

The most important general relativistic effect for X-ray line shapes is the gravitational redshift. The 
redshift formula is given in Eq. (22), and we intend to use that in our discussions of line profiles. How- 
ever, as remarked there, the formula is vaUd only in the case of stationary emitter and receiver, and the 
emitters in an accretion disk are not stationary. In fact, we expect relativistic velocities. In order to pre- 
serve the use of the simple equation for the redshift, we need first to know the photon energies as they 
would be measured by local stationary observers {dr/dT=d(j)/dT=d6/dT=Q), distributed in radius. We 
denote with the subscript Iso measurements made by such observers. Once we know the local photon 
energy distribution, it is then a simple matter to calculate the distribution at infinity by application of 
Eq. (22). Since disk material is, in the approximation that disk annuli describe circular orbits, passing 
the stationary observer with azimuthal velocity r d(j)/dtiso, we can simply apply a Lorentz transforma- 
tion that relates the spectrum in the disk frame to the spectrum measured by the stationary observer. 
The effective potential dictates a relationship between the angular momentum and the radius of circular 
orbits. Once we find this relationship, we proceed to the radial dependence of the velocity of circular 
orbits. 



5.6.1. Angular Momentum/Radius Relationship for Circular Orbits 

The condition for a circular orbit, dV/ dr = 0, yields a relationship between r and L. Using Eq. 
(29), this is 



m 



M 



m 



which has the two solutions 



r± 



{L/m) 
2M 



1± 1 



(L/m) 



1/2 



(30) 



(31) 



The stability condition d'^Vjdr^ > shows that a circular orbit at r+ is stable, while the one at r_ is 
unstable. Equation (31) thus provides the radii of both the unstable circular orbit and the stable circular 
orbit for a given value of L/m. We further see, that, owing to presence of the square root, the domain 
of allowed L/m yielding a circular orbit is restricted to L/m > 2\/3M. Setting L/m to this minimum 
value in Eq. (31), we find that the radius of the innermost stable circular orbit risco is equal to 6M. 

The existence of the innermost stable circular orbit, a uniquely general relativistic feature, can be 
appreciated by inspection of the right panel of Fig. 4, which shows the effective potential for several 
values of the angular momentum. The intersection of the dashed line with each effective potential curve 
marks the radial position of each stable circular orbit. The trend in the behavior of V with L/m is clear: 
as i/ro decreases, the extent of the potential barrier decreases, as well. As i/m — > 2\/3M, the barrier 
flattens, and disappears for L/m < 2\/3M, thereby eliminating the local minimum. All inward bound 
trajectories with L/m < 2^/2>M cross the event horizon. 

5.6.2. Energies of Circular Orbits 

From Eq. (28), setting dr/dr = for a circular orbit, 

2 



1 



2M 
r 



1 



{L/mY 



(32) 



From Eq. 30, we have 



m . 



r/M-2>' 



(33) 
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which, when substituted into Eq. (32), gives 

E _ 1 - 2M/r 
m ~ (1 - 3M/r)V2 



(34) 



This expression for energy is valid only for r > 6M because of the stability restriction. 

At the ISCO, where r/M = 6, we have {E/m)isco = (8/9)^/^. Therefore, the energy change per 
unit mass of a test particle as its circular orbit degrades from a large distance to the ISCO is 

A(':^')=l-f^V^'« 0.057 (35) 



mj \9 

i.e., the test particle loses an energy equivalent to about 6% of its rest mass energy as it works its way 
down to the ISCO in a succession of circular orbits. This is often taken as the maximum accretion 
efficiency of a Schwarzschild black hole, in the sense that iacc = V^accC^- 



5.6.3. Velocities of Circular Orbits 

We need to find the radial dependence of the velocities of circular orbits Vrj, = r d(f)/dt\so as they 
would be measured by locally stationary observers distributed in r. Proper time is measured in the disk 
frame, so that we can write 

# dr 

V0=r— — — , (36) 
dr dt\so 

and then relate dt\so and dr by a Lorentz transformation: 

dh,,^{l~vl)-^l^dT. (37) 
Using angular momentum conservation (Eq. 27), Eq. (36) becomes 

vl = ^ [L/mf (1 - vl). (38) 

Substituting for L/m from Eq. (30), and recalUng that that equation is valid only for stable circular 
orbits, some rearrangement gives 

which decreases monotonically with radius. The domain of validity is the same as the domain for which 
stable circular orbits can exist, i.e., r > 6M. The velocity of a particle in a circular orbit at the ISCO 
can be found from Eq. (39) by substituting r = 6M, from which we find that t;isco = 1/2. 

For the sake of interest, let us find the linear velocity and angular velocity f^oo of matter in a 
circular orbit as observed at infinity. To find the angular velocity, we start with the angular momentum 
conservation law (Eq. 27) and write 

^oo = ^ |, (40) 
at 

so that we need to relate dt and dr for a circular orbit. This is accompUshed by manipulating a simpU- 
fied form of the metric (Eq. 19), setting dr = 0, 

dr^ = (^1 _ dt' - dcj>^ (41) 
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which, after eliminating d^/dr, again using Eq. (27), becomes 



m] 



Substituting for L/m from Eq. (33) and simplifying gives 



"(ft 



1 - 



3M 



1/2 



(r > 6M). 



Substituting this back into Eq. (40), again using Eq. (27) gives 

1/2 

(r > 6M), 



MV 



r.3 ; 



which is identical to the Newtonian result. The hnear velocity is simply r^^, i.e.. 



1/2 



(r > 6M), 



which is also identical to the Newtonian result. 



(42) 



(43) 



(44) 



(45) 



5.7. The Motion of Light in the Schwarzschild Metric 

For Ught motion (technically, for the motion of massless particles), we derive the governing equa- 
tions by starting with the equations of motion for massive particles, then take the limit m — > 0. Starting 
with Eq. (28), we have 



1 + 



(Ljmf 



Using the equation for the conserved energy (Eq. 26), the previous expression gives 



~dt) 



2My 



2M 
r 



m 



'2 1 



,,2^2 



Taking the limit m ^ 0, and defining the impact parameter b = L/E,we find 



2My 



2M 
r 



For the azimuthal motion, we combine Eqs. (26) and (27) to give 



d(p E 
dt m 



\ r J 



which becomes 



1 - 



2M 



V 



r I 



(46) 



(47) 



(48) 



(49) 



(50) 



Therefore, once we know the impact parameter (see below), the Ught path can be calculated as pa- 
rameterized by our measure of time t. By eliminating dt between Eqs. (48) and (50), we can plot the 
trajectory r((^) according to 



dr 



1 



1/2 



(51) 
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Fig. 5. Effective potential for light motion near a Schwarzschild black hole (see Eq. 54). Three sets of (r, b) 
constitute the range of possibilties for light motion in the Schwarzschild potential, where b is the impact parameter. 
Light arrows indicate escape orbits, heavy arrows are capture orbits. The vertical dashed line shows the position 
of the unstable circular orbit. 



5.7.7. Effective Potential for Light Motion 

We wish to obtain an effective potential for light, analogous to that found for particle motion, in 
order to gain a quick qualitative understanding of light motion, without having to perform numerical 
integrations of the equations of motion. Following [152], we need the relations between time and radial 
displacements as measured by a local stationary observer and those measured by a distant observer, 
which results in 

Then, Eq. (48) can be written 



1 fdn^A^ _ 1 1 / 2M\ 
and we identify as the effective potential for light 



(53) 



^0 = 12(1- — )• (54) 



Note the simplicity of the description of light trajectories compared to particle trajectories, where, in 
the latter case, the geodesies depend both on energy and angular momentum. 

The effective potential for light is plotted in Fig. 5. The peak of V occurs at r = 3A/, for which 
Vpcak = 1/27M2. For the critical impact parameter, 6ciit = SVSM, driso/dtiso = at r = 3M. 
This is the only radius for which a massless particle can move on a circular orbit. However, the orbit is 
clearly unstable; any perturbation to the orbit will lead to either to capture or escape. 

For the purposes of determining the range of impact parameters that lead to capture by the black 
hole, the presence of the potential peak at r = 3M and the critical impact parameter divides space into 
two. First, consider light rays emanating from radii Vo > iM. From Fig. (5), we see that all outward 
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Fig. 6. {Left panel) Light trajectories in the Schwarzschild metric (see Eq. 51) for point of origin r = GM for 
twelve impact parameters confined to the range 5.0M-5.9M. {Right panel) The angle ^ between a light ray and 
a radius vector, as measured by a local stationary observer, used to define the impact parameter (see Eq. 56). 



directed photons escape to infinity. For inward directed rays, those for which b < 3^/3M are captured. 
For photons originating inside r — 3M, all inward directed photons are captured. Outward directed 
photons are captured if & > 3\/3M, and escape otherwise. 

Using the effective potential, we can calculate the fraction of photons that will be captured by the 
black hole if emitted isotropically by a stationary source at r. First, we need to know how to calculate 
b for various cases of interest. Suppose photons are being launched at various angles from a point r 
near the black hole. We can find b by referring to the azimuthal equation of motion. From Eq. (50), 
reckoning time according to a local stationary observer, we have 



dtiso 



V4, 



2M 
r 



1/2 



Using the coordinate system illustrated in the right panel of Fig. (6), this becomes 



r 1 



2AI 
r 



-1/2 



sin ip- 



Iso 



(55) 



(56) 



As an aside, a distinction needs to be made between angles measured by a locally stationary observer 
and those calculated by a distant observer. The relation between tpiso and the same angle according to 
a distant observer ip^o can be found as follows. Let the radial displacement measured by the LSO be 
denoted by ds. Then 



tani/jo 



dr 



r d(f> ds 
ds dr 



2M 



-1/2 



tani/j; 



1S0 5 



(57) 



where we have used Eq. (20) to relate ds and dr. 

Now, for r > 3M, we already know that all outward directed photons escape. Inward directed 
photons are captured if sinTAiso < 3V3{M/r)y^l — 2M/r. This results in 



/"capt c\ 



1- Wl-27( — 



1 - 



2M 



(r > 3M). 



(58) 
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Fig. 7. Capture fraction for stationary isotropic photon emitters plotted against radius for the Schwarzschild metric 
and for Euclidean space. In the latter case, the capture fraction is simply the fractional solid angle subtended at r 
by a sphere of radius 2M. For the Schwarzschild case, light bending is included. 



For r < 3M, we akeady know that inward directed photons are captured. Outward directed photons 
are captured if sin?/; > 3\/3{M/r)y^l — 2AI/r, which results in 



/capt — ^ 




(r < 3M). 



(59) 



The capture fraction is 1/2 for r = 3M, and approaches unity as r ^ 2M. A plot of /capt vs. r is 
plotted in Fig. (7). Calculating the captured fraction for material moving in a disk must also account for 
relativistic beaming. In other words, a moving isotropic photon source will, to varying degrees, appear 
anisotropic to the stationary observer, with its emission concentrated into the direction of motion. 



5. 7.2. The Speed of Light 

One of the more peculiar aspects of light propagation near a black hole is that those performing 
calculations in flat spacetime must account for the apparent reduction of the speed of light. According 
to a distant observer, the speed of light for pure radial motion, from Eq. (48), with 6 = 0, is 

dr , ^, 2Ar 
dt 



^ = ± 1- 



(60) 



which approaches unity for large r, and approaches zero as r 2M. For example, according to an 
observer at r, the time At for a light pulse to propagate radially from a position To < r to r is 



2 In 



r/2M - 1 



(61) 



M ' "' rol2M - 1 ' 

which shows the effect of the reduced speed of light (2nd term) compared to the Euclidean value ( 1 st 
term). Here we are measuring time in multiples of M. The conversion between seconds and centimeters 

is At(,) = {GM/c')M(,^y 
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Fig. 8. Time dependence of an initially radial planar light pulse near a Schwarzschild black hole. Time intervals 
are measured from initiation of pulse at t — 0. Deviation from perfectly circular wave fronts illustrates the 
combined effect of a variable speed of light and light bending, as calculated by a distant observer. Plots based 
on integrations of Eq. 5 1 . 



The apparent speed of light for pure azimuthal motion is, from Eq. (50), 



±- 



1 



± 1 



2M 

r 



1/2 



(62) 



where we have substituted for b from Eq. (56). Note that this is not the same as the speed of light for 
pure radial motion (Eq. 60). 

Combining Eqs. (60) and (62), we can evaluate the speed of light for arbitrary trajectories according 
to a distant observer 



dr\ ' 



dt J 



+ r 



dt 



1/2 



= 1- 



2M 



1 



2b^M 



1/2 



(63) 



While this apparent non-constancy of the speed of light might seem somewhat disorienting, it is im- 
portant to realize that we do not actually measure this speed. Measurements must be made locally. We 
can calculate what a locally stationary observer will measure by transforming Eqs. (60) and (62) to the 
frame of such an observer. 

1/2 



Iso 



dt 



Iso 





^ dr 
















'dt ^ 







dtu 



r ) 



dt 



1/2 



These two equations give 



1/2 



= 1. 



(64) 
(65) 

(66) 
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Therefore, locally stationary observers always measure c = 1. We, as distant observers, must contend 
with the fact that light appears to propagate at varying speeds, its speed depending on both the radial 
coordinate and the impact parameter. 

An example that illustrates the effect of a variable speed of Ught is shown in Fig. (8). A planar 
radial burst of light is emitted from r = 6M near a Schwarzschild black hole. Shown in the figure are 
the trajectories of a set of rays with a distribution of impact parameters, the endpoint of each trajectory 
indicating the progression of the leading edge of the pulse at four times: t = 6M, 12M, 18M, and 
24M, where the time is measured as a multiple of GM/c?. The reduced speed of light at small radii 
compared to larger radii is evident in the plot for t = 6M by the distortion from a perfect circle of the 
endpoints of the Ught ray pattern. With a knowledge of the disk structure (or, more likely, a model), one 
can imagine exploiting the unique light propagation physics in a black hole metric to locate the position 
of such a pulse as it irradiates different portions of the disk. Conversely, with a model of the origin for 
a pulse, or X-ray flare, the disk structure can be mapped, since the disk will respond to changes in 
the illumination pattern. Since such an effect occurs on a timescale that is measured in multiples of 
GM/c^, it is practical to observe "light echoes" in AGN (hundreds to thousands of seconds), but not 
in black hole X-ray binaries (~ 1 milUsecond). 

This leads to the concept of reverberation mapping, which is among the programs planned for 
future X-ray observatories, such as Constellation-X. We do not discuss the details of this ambitious goal 
here. However, a few comments are in order. Given that the model of X-ray fluorescence from AGN 
and black hole X-ray binaries involves the spectral response of an accretion disk to illumination by 
hard X rays, and that the hard X-ray continuum flux is observed to vary, it makes sense to ask whether 
changes in the shapes and intensities of fluorescence Unes in response to continuum variability can 
elucidate both the nature of continuum production and disk structure. The basics of such a program, 
for the case of a Schwarzschild black hole, have been laid out in [139]. It was shown that a variable 
source of localized hard X-ray illumination placed at the center of a disk would cause a characteristic 
response in the profile of a fluorescence line. In this case, with a centrally-located illumination source, 
the inner disk would respond first, affecting the red and blue extremes of the profile. Owing to light 
travel times, the rest of the profile would respond at later times. It is shown that, in principle, the mass 
of the black hole can be determined. Measurements of black hole spin (see §6) are also plausibly within 
reach [140]. In a more general case, reverberation could involve a non-axial illumination source in the 
Kerr geometry, time-dependent modifications to the disk structure, and multiple, overlapping X-ray 
flares, with a distribution in both intensity and position. In such models, the varying speed of light 
propagation figures prominently. For example, as discussed in [141], non-axial variable illumination 
causes the redward portion of the line profile to lag the rest of the profile, indicative of an inward 
moving pulse that sweeps across the disk, slowing to zero velocity at the horizon. 

Interestingly, current studies of AGN in which the X-ray continuum is observed to vary fail to show 
the expected variability in the iron Ka line complex [142] [143] [144], although efforts attempting to 
reconcile this problem are underway [145]. 



5.8. Minimum and Maximum Frequency Shifts 

Assume that a locally stationary observer measures the frequency distribution of a radiating particle 
that moves on a circular orbit at velocity and emits photons of frequency Vo in its rest frame. The 
minimum and maximum frequencies measured by the LSO are determined by a Lorentz transforma- 
tion: 

t'lso = li^o (1 ± v^,) (67) 
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where is given by Eq. (39), and 7 = (1 — Then from Eq. (22) we can find the frequency 

observed at infinity according to 



1/2 



This is straightforward, and works out to 



1 - 2M/r 
1 - 3M/r 



1/2 



(l-2M/r)i/2_^(M/r)i/2 



(68) 



(69) 



At the ISCO, we find i'^ = (\/2/3)i^o and — \f2vo [10]. This analysis also assumes that the 
observer at infinity sees the full red and blue shifts of the emitted photons. In more reahstic cases, the 
Une profile is narrower, since the disk is likely to be inclined. 

Note that, although an observer at infinity marks the velocity of a circular orbit as being = 
^jMjr, we cannot properly calculate the photon energies (as above) unless we analyze the problem in 
the local frame of the orbiting matter. The lesson is that there is no Lorentz transformation thai globally 
reduces curved spacetime to the Minkowski metric. 

We revisit the minimum and maximum frequency shifts in a later section, after we derive the 
analogous quantities for Kerr black holes. 



6. Spinning Blacl< Holes - The Kerr Metric 

We have already remarked on the ubiquity of angular momentum in the cosmos. Angular mo- 
mentum is also likely to play a role in the behavior of spacetime near black holes. Stellar precursors 
to Galactic black holes rotate, thus imbuing these stars with angular momentum. Near the center of a 
galaxy, gas that supphes fuel to supermassive black holes in AGN is in motion, and thus carries angular 
momentum with respect to the black hole. In both cases, then, we expect the black hole to possess an- 
gular momentum. Therefore, improved treatments of the physics in the black hole environment should 
include modifications to the spacetime geometry induced by spin [147]. 

In 1963, Roy Kerr [148] presented his solution for the metric outside a spinning object at the First 
Texas Symposium on Relativistic Astrophysics (see [149] for a description of these proceedings). The 
physical imphcations of this solution were pursued over the course of the next few years. More recently, 
observable properties of the emission and timing properties of accretion disks have led to attempts to 
measure black hole spin. 

The Kerr metric is given in Boyer-Lindquist coordinates [150] as 

dt^ + ^^"7'"'^ dtd<^-\ dr' -Y.d6^- Rle sin' d<i>\ (70) 

where 

Y, = r^ + a? cos^ 9 




A = - 2Mr + 0? (71) 

Ri g = + 0,2 + 2Mra^Tr^ sin^ 9. 
For motion in the equatorial plane {6 = 7r/2), the Kerr metric is 

dr^=(l-'^\ df + ^M^dtdcl>-Cdr^- Rl #^ (72) 
\ r J r A 
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where 

Rl=r^ + + ——_ (73) 
We show below that Ra is the reduced circumference in the equatorial plane of a Kerr black hole. 

6.1. The Kerr Event Horizon 

The radial coordinate of the horizon rjj is defined as that radius at which the coefficient of the dr"^ 
term of the metric blows up. For the Kerr metric this occurs when A = 0, which gives 



th = M + - a^. (74) 

This shows us that the maximum spin consistent with a real- valued radial coordinate for the horizon is 
simply M. The radial coordinate of the horizon thus ranges from 2M (o = 0) down to M (o = M). 
Detailed considerations concerning the spin-up of a black hole by accreting matter shows that the 
maximum spin is not M, but rather 0.998M [151]. The case a = M (or a = 0.998M) is referred to as 
an extreme Kerr black hole, and the black hole is said to be maximally spinning. 

6.2. The Static Limit and the Dragging of Inertial Frames 

The following derivation follows that in [10]. The relation gnuP^p" = —w? is equivalent to the 
condition Qfj^^u'^u" = — 1, where u** is a component of the velocity four- vector. Consider an observer 
at fixed r and 6, so that dr = and d6 = 0. Then 

gtt {u'f + 2gt4, u'u'^ + g^^ {u'^f = -1, (75) 

where the metric coefficients can be read off from Eq. 70. Rearrangement of the previous equation 
gives 



t\2 



K) 



9tt + '^gt4> —r + 



= -!• (76) 



If we let m"^/ m* — d(j)/dt = co, then we require 

^2^2^a; + ^<0. (77) 
94>4> 94"t> 

This restricts ui to the range 

1 5*0 1 



1 



9ttg4>4, 
9t<t> . 



9<t><t> 




9tt94,<t> 



9U 



(78) 



or Wmin < oj < Wmax- Notc that is always negative, and that gr^^ is always positive. The metric 
coefficient gu can be positive or negative: 



gtt < if r>M + y/M'^ - a'^ cos^ 6 

(79) 

gtt>0 if M <r < M + VM^ - cos^ Q. 

Thus whenever gtt > (small radii), we see that Wmin > 0; the observer is swept around the black 
hole, and no action on the part of the observer can change that. Where gu < (large radii), Ui-ain < 0, 
and retrograde motion is allowed. The critical boundary, the mathematical surface inside which matter 
is dragged irresistibly in the prograde direction is called the static limit, which, from Eq. 79 is 



rstat = M+ VM2-a2cos2 e, (80) 



© NRC Canada 



32 



unknown Vol. , 



spin vector 
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plunging 
region 
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Fig. 9. Schematic of a disk-accreting Kerr black hole (a — 0.9A/) in cross-section, using t, r, 9, cj> coordinates, 
illustrating various geometric relations. The dark circle is the black hole, bounded by its event horizon. The gray 
region, delineated by the horizon and the static limit {dashed lines) is the ergosphere, cylidrically symmetric with 
respect to the spin vector. The accretion disk is assumed to lie in the equatorial plane of the hole. The disk is 
truncated near the hole according to the radius of the innermost stable circular orbit appropriate to the spin. Inside 
the ISCO, matter free falls through the plunging region and crosses the horizon. Tapering of the disk is schematic 
only, suggesting a general trend of geometrical thickening with radius. Also, see Fig. (12). 



which corresponds to gu — 0. 

The peculiar behavior of spacetime near the static limit can be illustrated [152] if we imagine 
launching light pulses in the azimuthal direction, first in the prograde direction, then in the retrograde 
direction, at some radius r near a spinning black hole. We wish to find d(f)/dt, the angular velocity 
of the light pulse according to distant observers, for these two cases. For simplicity, we work in the 
equatorial plane. Then, with 6 — tt/2, dr = 0, dO = 0, and the condition Qf^vU^u^ = 0, we find the 
following quadratic for d(j)/dt, 



From Eq. 80, the static limit in the equatorial plane is Tgtat = 2Af , which, when substituted into Eq. 
82, gives the two solutions d<l)/dt — AMa/rR^ = a/{a'^ + 2AI^) and d4>/dt — 0. The pulse emitted in 
the retrograde direction appears to stand still. Inside the static limit, retrograde pulses appear, in fact, to 
move initially in the prograde direction. Although we do not pursue it here, it is worth noting that the 
region between the static limit and the horizon is called the ergosphere, after the Greek word ergon, for 
work, since energy can, in principle, be extracted from it [153]. An illustration of the size and shape of 
the ergosphere is shown in Fig. (9), along with other components of an accreting black hole. 

An expression for the rate at which spacetime is dragged around a spinning black hole can be 
derived from the momentum conservation law, which, as in §6 for the Schwarzschild metric, we derive 
from the Euler-Lagrange equation for cf). A suitable Lagrangian for motion in the equatorial plane is 
given by 




(81) 



which has the solution 




(82) 




(83) 
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Fig. 10. Plots of photon trajectories with zero angular momentum near a Kerr black hole, for a — 0.2M (left 
panel) and a = 0.998M (right panel) illustrating the frame dragging effect. Circles are lines of constant r. Plots 
are based on integrations of Eq. 116. 



where the dot denotes differentiation with respect to any parameter p along the particle's world line 
[10]. The Euler-Lagrange equation for the (/)-coordinate results in 



(84) 

dp r dp 

If dp — dr/m, then we find an expression that can be compared to the case for the Schwarzschild 
metric (Eq. 26), 



2aM dt 



Rl^^- + ^^- (85) 
dr m r dT 

to which it reduces for a 0. The second term on the right, which is proportional to the black hole 
spin, shows that a particle with d(l)/dT = must have negative angular momentum. In fact, for massive 
particles and massless particles alike, if L = 0, then 

d0 2aM _ 

d^^7i?r="'^^^«- ^^^^ 

According to a distant observer, a particle with zero angular momentum appears to rotate in the direc- 
tion of the black hole spin (i.e., following the right-hand rule). This effect is known as the dragging 
of inertial frames. The drag frequency increases with decreasing radius. Frame dragging can also be 
described as follows: Imagine yourself in a rocket near a spinning black hole. Looking at the fixed stars 
for reference, you find that you require a thrust in the positive azimuthal direction (i.e., pointed toward 
the spin direction) to hold the rocket in a position such that the stars overhead appear not to move. It is 
said that spacetime itself is moving. 

The dragging of inertial frames occurs outside of any spinning body. As an example, we use Eq. 
(86) to estimate Wdrag caused by the Earth's rotation. In conventional units, the drag frequency at the 
Earth's surface is uidiag ~ 2GIuje / {<? R%), where / is Earth's moment of inertia. Re is its radius, 
and loe is its rotation frequency. The moment of inertia of the Earth is Q.Z'iMER% [8], where the 
deviation from 2/5 in the multiplier is due to the non-uniform density distribution of Earth's mass. 
Taking Me = 6.0 x lO^"^ g. Re = 6.4 x 10*^ cm, and uje = 7.3 x lO'^ s-\ we find Wdrag = 219 mas 
yr~^, large enough to be measured with current technology. Recently, the LAGEOS mission, a pair of 



L 2aM dt 
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laser-ranged satellites on orbits with aphelia of approximately 12,000 km, measured Wdrag> finding a 
value that is 99 ± 5 percent of the value predicted by GR [154]. Gravity Probe-B is expected to reduce 
the uncertainty in the measurement to 1%. 

Interestingly, we find for the Earth that a = 7.3 x 10^ Me- The restriction a < M is clearly not 
a general constraint on material bodies; instead, it can be viewed as a necessary condition for matter 
lying inside its event horizon. 



6.3. Gravitational Redshift in the Kerr Metric 

In the Schwarzschild metric, we found the relation between proper time and coordinate time by 
considering two stationary observers at different radii, and found that the tt-component of the metric 
contained the necessary information. An analogous approach in the Kerr metric involves an observer, 
called a zero angular momentum observer (ZAMO), moving along the 0-direction with an angular 
velocity Wdrag given by Eq. 86. By a simple transformation, the metric can be diagonaUzed, i.e., we 
can decouple the t and cf) coordinates, then relate dt to dr. Let d(t) = Wdrag dt + d(pr, where d(pr is 
an azimuthal displacement measured in the ZAMO frame. The case of greatest interest is that of an 
emitter in the equatorial plane. Therefore, for pure azimuthal motion, Eq. (72), with dr = 0, becomes 



dr' 



{l-'-f-)df + '-^dtd^-Rld^\ (87) 



Substituting the d^ transformation into this form of the metric yields a new metric that uncouples the 
t- and ^-coordinates: 

dr^ = (^1 - — + -^jj^j dt" - Rl dct>l. (88) 

This shows that the reduced circumference in the equatorial plane of a Kerr hole, the analog of r used 
to describe the Schwarzschild geometry, is Ra- Equation (88) gives us the relation between time for a 
ZAMO (set d(j)r = 0) and a distant observer, which we write as a redshift formula for frequency: 

2M 4MV\'/^ 



or, after expanding Ra, and rearranging, 

(90) 



r'^+a'^ + 2Ma?/r 

which reduces to Eq. (22) in the limit a ^ 0, the Schwarzschild case. As a consistency check, note 
that the surface of infinite redshift, the outer solution of -\- a?' — 2Mr = 0, is identical to our earlier 
result, Eq. (74). 

The multiplier of i^o in Eq. (90) is sometimes called the lapse function. Since we use this function 
later in the calculation of orbital velocities, it is defined here and given the label a: 



\/r^ +0? - 2Mr 

a = . (91) 

6.4. Circular Orbits in the Kerr Potential 

Evaluation of the Euler-Lagrange equation for i provides an energy constant of the motion, 

/ 2M\ . 2Ma ■ ^ 

1 ]t+ = (92) 

\ r J r 
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Fig. 11. (Left panel) Radius of innermost stable circular orbit plotted against the spin parameter a*. Radii of the 
ISCOs range from QM (a, = 0) to M (o* = 1). For physically realizable Kerr black holes, the maximum spin 
is a, = 0.998, corresponding to an ISCO radius of 1.2AM. {Right panel) Dimensionless accretion efficiency rj 
plotted against a*. Efficiency ranges from 0.057 (a* = 0) to 0.42 (a* = 1). For a, = 0.998, 77 = 0.34. 



using the Lagrangian in Eq. (72). We have akeady evaluated the <j) component (Eq. 84), which we 
reproduce here, 

+ = (93) 

r 

These two equations can be solved for t and (f) in terms of the energy and angular momentum, yielding 

. ^ ErRl 2LMa ^^^^ 
rA 

. L{r - 2M) + 2EMa 

0= 7 • (95) 

rA 

Substitiution of Eqs. (94) and (95) into the Lagrangian (Eq. 83), noting that A — —m? 12, yields an 
equation for r, which we write in terms of dr / dr, substituting p ^ r/m, 

r^(±\\s^(^\ \rI - AMa f A V^) _ _ 2M) ~ rA. (96) 

\dT J V"^/ \m J \m J \"^/ 

Stable circular orbits can be found by imposing three conditions on S": (1) S* = (turning point); 
(2) dS/dr = (circle); (3) d^S/dr^ < (stability). The first two conditions yield a system of two 
equations that can be solved for E/m and L/m, restricting those quantities for a given r and a: 



E r - 2M + aJM r 

— ^ , , (97) 

m (r2 - 3Mr + 2aVMr) 1/2 

L_fMV'^ r^-2aVW + a^ 
m^yrj (r2 -3Mr + 2a\/Mr)i/2' 

The third condition requires that 

eV . 2M 



- -l + ^<0. (99) 
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Substituting from Eq. 97, and defining the dimensionless spin parameter a* = a/M, the previous 
restriction becomes 

X 2 / X 1/2 

tt) -S^T + Sa*!^) -3a2 >0, (100) 

where the equality can be treated as a quartic equation in {r/MY/"^. The quartic has been solved [146] 
for the radius of the innermost stable circular orbit for a specified spin parameter: 



' ISCO 



= M 



3 + ^2 - V(3-^i)(3 + ^1+2^2) 



(101) 



where Zi and Z2 are defined by 
Zi = 1 + (1 - a2)V3 [(1 + a*)V3 + (1 - a,)V3] 

(102) 

Z2 = (Z? + 3a2)V2. 

Equation (101) generalizes the result based on Eq. (31) for the Schwarzschild metric. As shown in he 
left panel of Fig. (11), rigco is a monotonically decreasing function of spin, ranging from a maximum 
of 6M for a Schwarzschild hole to a minimum of M for a maximally- spinning Kerr hole. In Fig. (12) 
we show schematically the variation of the position of the ISCO of an accretion disk compared to the 
size of the event horizon and the shape of the ergosphere, where, for example, it is seen that for large 
values of a* the ISCO lies inside the ergosphere, and approaches the event horizon. 

If we assume that the accretion efficiency 77 is determined by the energy of the innermost stable 
circular orbit, then the maximum efficiency occurs for the case of maximal spin. With r = M and 
a = M, Eq, (97) gives E/m = 1/a/3. Therefore, for a disk extending down to nsco for a maximally 
spiiming black hole, the accretion efficiency is 

via. = 1) = a(^^^ =l-(^iy « 0.423, (103) 

which can be compared to Eq. (35) for the Schwarzschild metric. The energy conversion efficiency t] 
is plotted as a function of the spin parameter a* in the right panel of Fig. (11). 

The velocity of matter in a circular orbit around a Kerr black hole, as measured by a ZAMO, can 
be expressed in terms of the orbit's angular frequency as reckoned by a distant observer (O) as 

R 

V4, = {fl - Wdrag), (104) 

where Ra (Eq. 73) is the reduced circumference of a spinning hole, Wdrag (Eq. 86) corrects for the 
dragging of inertial frames, and the lapse function a (Eq. 91) accounts for the transformation between 
measures of time. An expression for O [155] [156] is most easily found by starting with the equation 

of motion 

where the Christojfel symbols are given in terms of the metric g^v and the inverse metric 5'"' by 



dx^^ dx^ dx'^ 



(106) 



whose derivation can be found in [137]. For a circular orbit, using Eqs. (94) and (95), the first term on 
the left of Eq. (105) vanishes for all A. If we choose A = r, then 

dx'^dx" = Tl, dt^ + 2 r;, dt dcP + d<P^ = 0, (107) 
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Fig. 12. Comparison of the radius of the event horizon (solid circles), the static limit {dashed lines) and ergo- 
sphere, and position of the innermost stable circular orbit for black holes of various spin parameters. 



or, letting ft — dcf)/ dt, 

r;^f]2 + 2r;,f] + r;, = o. (i08) 

An evaluation of the Christoffel symbols gives 

r[, = ig-(2A//r2) 

= i5-(-2Ma/r2) (109) 

= ig- (2MaVr2 - 2r), 
which, when substituted into Eq. (108), leaves us with the simple quadratic equation 

Ma^ \ 2Ma „ M 



r]n^-^^n+'—=0. (110) 
The solution corresponding to prograde motion is 



jM/r 1 1 
O — X ' — dill 

r + aVMA M (r/A/)3/2 + a, ■ ^ ^ 

This expression for can also be derived by taking the ratio of and i from Eqs. (95) and (94), 
respectively, then substituting for E and L from Eqs. (97) and (98), respectively, but the algebra is 
tedious. 
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Fig. 13. {Left panel) Azimuthal velocities for equatorial circular orbits as a function of radius for five values 
of black hole spin, as measured by locally stationary observers. Curves are labeled by the spin as fractions of 
M. Vertical lines indicate radii corresponding to the innermost stable circular orbit. (Right panel) Azimuthal 
velocities for a selected group of spin parameters (as labeled), illustrating the development of a local minimum for 
a. > 0.9953. 



Substituting Eqs. (73), (86), (91), and (111) into Eq. (104) yields an expression for the velocity of 
circular orbits [146], 



which reduces to the Schwarzschild case, Eq. (39), for a = 0. For the other limit, a — s- M, the 
numerator can be factored to give 

MV''^ ^3/2 ^ ^^1/2^ _^ ^,^^1/2 „ ^^3/2 



\r) (r + M3/2r-i/2)(ri/2+Mi/2)' 

which shows that as r ^ M, the velocity approaches — 1/2, which is, perhaps surprisingly, 
identical to the velocity at the ISCO of a Schwarzschild black hole. The left panel of Fig. (13) shows 
the radial profile of i'0(r) for several different spin parameters. The maximum orbital velocity attained 
is v'^'''' = 0.624 for a, = 0.9268, which occurs at the ISCO with n^co = 2.14M. 

Another feature worth noting is that for the extreme a = M case, the maximum velocity is ap- 
proached from above, i.e., v^{r) is not monotonic. For more realistic cases, where the maximum value 
of a, is 0.998, the velocity at the ISCO is always approached from below. But this velocity is not 
always the maximum velocity. For = 0.9964 and above, the radius of maximum velocity r('i;niax) 
is distinct from the ISCO, as can be seen in the right panel of Fig. (13) (lowest two curves), with the 
largest separation r (wmax) — ''isco occurring for the maximally spinning case, where r(timax) = 1 .72M 
and where Wmax — 0.573. The non-monotonic behavior of (r) sets in for a* > 0.9953. For a* in this 
range, ^^(r') has a local minimum that moves inward as the spin parameter increases. Recently, it was 
proposed [157] that this feature of may be related to the fixed ratio of high frequency quasi-periodic 
oscillations (QPOs) observed in Galactic black hole X-ray binaries. 

Having obtained formulae for the gravitational redshift (Eq. 90) and the velocities of circular orbits 
(Eq. 112), we now follow the same procedure used for Schwarzschild black holes (§5.8) to find the 
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Fig. 14. Minimum and maximum normalized line energies observed at infinity vs. disk radius for several black 
hole spin parameters, assuming circular orbits and 90 degree inclination. The upper set of curves shows the maxi- 
mum blueshifts for the eleven spin parameters a, = [0., 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.998], working 
downward from the dashed line (Schwarzschild case). The lower set of curves shows the corresponding maximum 
redshifts, working upward from the a, — case. Pairs of curves delineate the maximum width of the line pro- 
file for a given radius, neglecting thermal Doppler widths, and any possible radiative transfer effects. Curves are 
truncated at the ISCO corresponding to each spin parameter. 



minimum and maximum line energies as observed at infinity as a function of radius for several spin 
parameters. The result is shown in Fig. (14), where the extrema are plotted against radius for eleven 
values of the spin parameter From that figure one sees that the differences between spinning black 
holes and a Schwarzschild black hole are fairly minor for r > 6M, viz., the maximum extent of the 
blue wing is the same to within about 10%, for r > 6M, and the maximum redshift is even less 
distinct when comparing different spins. If, for example, the disk of a maximally spinning hole did 
not, for some reason, produce line emission inside of about 6M, one would be hard-pressed to discern 
spin by measuring the line width alone. Figure (14) also shows that the line width is bounded by 
the Schwarzschild "envelope" for r > 6M (dashed line). The maximum blueshift is a monotonically 
decreasing function of spin; the maximum blueshift never exceeds the spin-zero case E/Eq = 1.41 at 
r = 6M, which may be surprising. For a given spin, and moving to smaller radii, the modest change 
in the maximum blueshift illustrates the trade-off between increasing orbital velocity and increasing 
gravitational redshift. For large spin values, the maximum blueshift is actually redshifted from line 
center, showing the dominant efffect of gravitational redshift. Thus the radius of maximum blueshift 
is not always the innermost stable circular orbit. This is most obvious for maximally spinning holes, 
where the energy of the blue side of the line approaches w 0.2Eo, since the ISCO approaches the event 
horizon. The maximum redshift always occurs at the ISCO, with the minimum value of E/Eo being 
0.060 for a* = 0.998. The maximum redshift increases monotonically with spin, since the effects of 
increasing velocity and increasing gravitational redshift reinforce each other, rather than offset each 
other, as they do when evaluating the maximum blueshift. One also sees from the figure that the total 
line width changes with radius. For example, an emission profile from an annulus near r — 2M in 
a maximally spinning hole would appear relatively narrow, with a highly redshifted centroid energy. 
The narrowing of the line follows from the overwhelming effect of gravitational redshift as the ISCO 
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approaches the event horizon. 

In any realistic situation, one expects fluorescence from a large range of radii. A model hne profile 
is built up annulus by annulus, where account is taken of the intrinsic surface brightness variation, 
the fraction of emitted photons that escape the system and find their way to the observer (given light 
bending), and the dynamical aspects discussed above. These model profiles can then be compared 
to spectroscopic data, which, in principle, allows one to extract several parameters of interest (disk 
inclination, radial emissivity profile, disk inner edge, etc.). Considered as a diagnostic of black hole 
spin, the distinguishing characteristic of line emission is, from Fig. (14), the extent of the red wing. 
By contrast, the blue wing is more sensitive to the inclination than to the spin. Given an observed Une 
profile, the inclination, to a first approximation, can be determined by the position of the blue edge. 
The spin can then be estimated from the redward extent of the hne profile. It's not quite that easy, 
unfortunately, but computer programs allowing fits to relativistic line profiles exist and are in wide use 
(e.g., [21] [158]), made available, for example, as part of the XSPEC [159] or ISIS [160] data analysis 
packages. 

The detection of highly redshifted iron Ka emission has led to inferences of non-zero spin, and 
has given rise to the idea that X-ray lines can be used to measure black hole spin. A dramatic example 
of a red-shifted iron Ka hne profile has been documented for the Seyfert 1 galaxy MCG-6-30-15 
[161] [162], in which the line profile was observed to vary dramatically, from a profile for which the 
emission appeared to be dominated by large radii, and then to a profile in which the bulk of the emission 
was observed to lie below the rest energy. This latter result implies that the emission originated from 
r < 6M, which, given the assumption that radiation from within the ISCO is unobservable, imphes 
black hole spin. In fact, a fit to the data requires a* > 0.94 [163]. 



6.5. The Plunging Region 

As orbiting matter in accretion disk reaches the ISCO, it is typically assumed to lose all rotational 
support and free fall from there through the event horizon. In terms of test particle trajectories, the 
gradually shrinking circular orbits are replaced by a rapid spiraling infall. We can calculate such a 
trajectory by eliminating dr between Eqs. (27) and (28), from which we find the following expression 
for the trajectory of a test particle in the Schwarzschild metric: 

acp L/m l\m J \ ^ J\ /. 

If we set E/m = a/8/9 and L/m = 2v^M, the values for a circular orbit at the ISCO, then add a 
slight radial displacement inward, an orbit such as that shown in Fig. (15) results. 

The sudden increase in radial velocity results in a reduced disk optical depth, which may reduce 
the X-ray albedo. The free falling matter is also conventionally assumed to be stress-free, so that a hot 
corona cannot be supported within the ISCO. Therefore, irradiation inside the ISCO by hard X rays 
is reduced relative to matter at larger radii. Finally, the drop in density may imply an increase in ^, 
which reduces the overall efficiency of fluorescence hne production. Taken together, these aspects of 
disks imply that any fluorescence must occur outside the ISCO. In other words, the radii over which 
observable fluorescence emission originates has a hard iimer limit of r = 6M for a Schwarzschild black 
hole. Given such a constraint, the assumption of a zero-spin black hole therefore Umits the maximum 
redshift of the hne profile (see Fig. 14). Profiles violating this limit have thus been taken to imply spin 
(see §6.4). 

The whole concept of a clearly demarcated iimer disk edge has recently been called into question, 
based both upon test particle trajectories [164] and upon three-dimensional MHD simulations of ac- 
cretion disks [165]. In the former paper, with a source of hard X-ray illumination located on the axis 
above the disk plane, it was found that substantial fluoresecence hne flux from inside the ISCO may 
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Fig. 15. Trajectory of free falling test particle inside the innermost stable circular orbit — the plunging region 
— of a Schwarzschild black hole (see Eq. 114) according to an observer at infinity. The energy and angular 
momentum are set to the corresponding values at r = 6Af . In order to perturb the particle off the stable orbit at 
r — 6M, the initial radius was set to 5.99M. 



emerge from the system, and that this component of the flux appears at redshifts exceeding the usual 
limit imposed by a hard cutoff at r = 6M, thus mimicking emission from a spinning black hole. The 
importance of emission from inside the ISCO of a non-spinning hole has been contested [140], how- 
ever, based upon the absence of an iron K absorption edge in the well-studied source MCG-6-30-15 
that is predicted by modeling calculations. In [165] stress is found to be continuous across the ISCO, 
invalidating the usual assumption that the stress vanishes there. Moreover, it was found that the matter 
density distribution need not conform to the simpler model wherein the density decreases monotoni- 
cally with decreasing radius, but that dumpiness may permit the survival of matter with lower values 
of ^. Thus there is some doubt cast upon our ability to make clean measurements of black hole spin 
using the iron Ka profile alone. 

6.6. The Motion of Liglit in tlie Kerr Metric 

A general treatment of light motion in the Kerr metric is considerably more complex than the 
Schwarzschild case (see, e.g. [166]. For example, except for the case of motion in the equatorial 
plane, trajectories are non-planar (see Fig. 16). Discussions of the techniques used to calculate photon 
geodesies near Kerr black holes may be found in [166] [167] [168] [169]. 

We can, however, work out the relatively simple example of a massless particle with L — moving 
in the equatorial plane of a spinning black hole. We find an expression for dr/dcj) from the following 
two equations: 



(1 - 2M/r) dt^ + (4Ma/r) dt dcj) - (rVA) dr^ - Rl dcj)"^ = 0. 
The first of this set comes from Eq. (84) with L = 0, and the second is simply the metric (Eq. 72) with 




d(t)- {2Ma/r)dt = 



(115) 
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Fig. 16. Trajectory for a massless test particle near a maximally spinning black hole, showing the generally non- 
planar character of photon trajectories in the Kerr spacetime. The path is projected onto the x-z plane, in which 
the angular momentum vector of the black hole lies. Photon is launched from r — 1.27-Rg, with an initial mo- 
mentum vector predominantly in the positive y-direction, with a small positive x component and a small negative 
z component. The photon executes several "circuits" before impacting the horizon near the pole. The calculation 
was performed following the methods described in [166]. 



dr = 0. Equation (115) is a system of two equations in the variables dt, dr, and d(f). By eliminating dt, 
we arrive at an expression involving only dr and dcf) that describes the photon trajectory: 



dr ^ 



2Ma' 



(116) 



for which a straightforward numerical integration produces the plots shown in Fig. (10), where trajec- 
tories for L = photons follow the spacetime drag for two values of the spin parameter 



7. X-ray Fluorescence Spectroscopy of Accreting Black Holes 

Having introduced the basic aspects of GR and accretion disk theory, we look now at the production 
of X-ray fluorescence line emission in black hole accretion disks. We start with a brief discussion of 
photoionization codes, and introduce the ionization parameter concept. Then we treat the phenomenon 
of "reflection," the spectral response of an optically thick medium to irradiation by a hard X-ray source. 
We close the section by presenting a few aspects of fluorescence at the atomic scale. The treatment of 
the latter topic is primarily from the point of view of atomic modeling, rather than from a quantum 
mechanical standpoint. 



7.1. X-Ray Photoionization Codes 

Photoionization codes, such as those described in references [172] and [173], are used to determine 
the effect of a radiation field on a gas of specified chemical composition, and the self-consistent effect 
that passage through the gas has on the radiation field. In other words, the opacity determines the effect 
of the gas on the radiation field, but the radiation field partly determines the opacity, primarily through 
its influence on the charge state distribution (ionization balance) [174]. In addition to acting as the dom- 
inant agent of ionization, photoionization also heats the plasma, since suprathermal photoelectrons are 
thermalized after interacting with the local population of MaxweUian-distributed electrons. Compton 
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scattering and the Auger effect (see below) also contribute to plasma heating. One may add a source 
of non-ionizing heating to mimic the sum of various other processes [175]. Plasma cooling is primar- 
ily through recombination, coUisionally-excited line emission, bremsstrahlung, and inverse Compton 
scattering. Also, imposing the constraint of charge neutrality controls the overall free electron to ion 
ratio. The explicit calculation of heating and cooling thus couples the energy equation (radiative heat- 
ing = radiative cooling) to the ionization equations and the neutrality equation, constituting a complex 
system of equations, whose solution requires an iteration scheme. 

Three approaches to calculating level populations can be used: (1) the nebular approximation, in 
which only the ground state has a significant population, (2) the Saha-Boltzmann approach, in which 
the populations of excited levels k are given by = {gk/gi) cxp{—Ek/kTe), where g denotes the 
statistical weight factor, Eh is the level energy with respect to ground, and kT^ is the local electron 
temperature, and (3) detailed level accounting, where level populations are calculated explicitly by 
diagonalizing the rate matrix, the elements of which include all rates into and out of each level, thus 
requiring many thousands of atomic rate coefficients. Once the level populations are specified, the 
local contribution to the overall spectrum is determined. Finally, by more or less approximate methods, 
radiation transport of lines and continua from their sites of creation to the observer is accounted for 
[176] [177] [178] [179]. Thorough discussions of photoionization codes and their apphcations can be 
found in [172J L180j |181J [182] [183]. 



7.7.7. The Ionization Parameter 

The physical conditions in X-ray photoionized plasmas are, for a given ionizing spectrum, often 
described in terms the ionization parameter [184]. This quantity arises naturally from the steady-state 
equations of ionization equilibrium. Let /3j, Ci, and Ofi+i denote the photoionization rate (s~^) of 
charge state z, the collisional ionization rate coefficent (cm-^ s~^) of i, and the recombination rate 
coefficient (cm^ s~^) of charge state i + I, respectively. The term a^+i accounts for all two-body 
recombination processes. Ignoring three-body recombination, as well as coupling to charge states more 
than one charge away, the steady state equations of ionization equilibrium can be written 

[fii +neCj(Te)] Ui = Heai+iiTe) n^+i, (117) 

where the rate coefficients for recombination and collisional ionization depend explicitly on the elec- 
tron temperature Te- In terms of the photoionization cross-section ai and ionization threshold energy 
Xi of charge state i, the photoionization rate for a point source of ionizing continuum can be written as 
an integral over photon energy: 

L '^'^ 



A = ^ / dE^a.iE), (118) 



Xi 



where Se^s, the spectral shape function, normahzed on a suitable energy interval. Denoting the integral 
in Eq. (1 18) by Eq. (1 17) becomes 

m + Ci{Te) Hi = ai+i(Te) rii+i. (119) 



Let ^ = Lx/nev'^, which is called the ionization parameter. Then 



(120) 



Ui Q!i+i(Te) Q!i+i(Te)' 

where the second approximate equahty apphes for many cases of interest. Such a plasma, where the 
ionization and energetics are dominated by the influence of an X-ray field, is referred to as an X-ray 
photoionized plasma. More generally, when the source of X-ray illumination is not a point source, 

l = i^^ .12,) 
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Fig. 17. Iron charge state distribution vs. ionization parameter. Figure from Kallman et al. [185]. 



where F is the energy-integrated, angle-averaged, ionizing flux, and n is the particle number den- 
sity. There are other, related forms of the ionization parameter in use (e.g., [175]), but we will use ^ 
throughout. In terms of ^, the photoionization rate corresponding to a single ionic cross-section is 

A = < dE a,{E). (122) 

For the X-ray sources of interest, typical values of ^ range from approximately 1-10'' erg cm s^^. 
Since the energy equation and the equations of ionization balance are coupled, the output from the 
calculation of the physical state of a gas using a photoionization code gives both the temperature and 
the charge state distribution as a function of ^. A pertinent example is given in Fig. (17), which shows 
the distribution in ^ of all 27 charge states of iron [ 1 85] . 



7.2. X-ray Fluorescence Lines 

As discussed above, in accretion-powered X-ray sources, such as AGN and X-ray binaries, repro- 
cessing of a hard X-ray continuum in relatively cool matter (10^ - 10® K) can generate intense iron K 
radiation from more neutral iron species of relatively low charge [83]. Iron fluorescence can be quite 
prominent, and, in fact, has long been known to constitute an essential component of disk spectra [186] 
[187]. 

The reprocessing mechanism, for any species with more than two bound electrons, begins with 
photoionization of a Is electron by a photon with energy e above the K edge, sending an element A in 
charge state z — 1 to charge state i, where i is in a quasi-bound state that is coupled to the continuum, 
denoted by the double asterisk below: 

e + Ai-i A** + (K shell photoionization). (123) 
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Fig. 18. Schematic of atomic processes involved in production of photoionization-driven K spectra. Three charge 
states, i — and i + 1 are represented. Level 1 is the ground state of i, level 2 is an excited level, levels 3 and 4 
are autoionizing, lying above the first ionization limit of i (dashed line). The low-lying level b in i — 1 is populated 
by photoexcitation from the ground level a, and decays radiatively back to a (wavy lines). (In the more general 
case, coUisional excitation and de-excitation also affect the population distribution.) K-shell photoionization from 
levels a and b populates levels 3 and 4, respectively, which can decay radiatively by Ka emission to levels 1 and 
2, respectively. Alternatively, levels 3 and 4 can autoionize, leaving as products an Auger electron and ion i + 
The relative intensities of 3 ^ 1 and 4^2 thus depend on the level population distribution in 1 — i. The escape 
probabilities of these two lines depend on the level population distribution in i. 



The intermediate state A**, since it lies above the ionization threshold, has a non-vanishing e^- 
Coulomb repulsion matrix element with a final product state, part of which consists of a continuum 
electron. Therefore, there is a non-vanishing probability assigned to the reaction 

A** ^ Ai^i + (autoionization). (124) 

The ejection of an electron by this mechanism is called autoionization. The configuration or energy 
level associated with A** is called an autoionizing configuration or autoionizing level, respectively. 
Referred to as the Auger effect for Is vacancy states, autoionization is often the dominant decay route 
for the Is-hole state. The reaction products are an Auger electron, i.e., an electron with a kinetic energy 
that is characteristic of the atomic energy level structure, and an ion in charge state i + 1, where, in this 
example, we assume for simplicity that i + lis left in its ground state. Still, it is worth mentioning that 
autoionization does not always leave a product ion in the ground state. It may be excited, and it can even 
be autoionizing. In the latter case, several charge states are coupled by autoionization; a single K-shell 
photoionization thus leads to a vacancy cascade, and must be accounted for when calculating the charge 
state distribution [188]. An extreme case relevant to astrophysics occurs for K-shell photoionization of 
neutral iron (Fe I), which initiates a vacancy cascade that can leave Fe X [189]. However, even in the 
simplest case described above (Ai^i Ai — > ^i+i), three charge states are coupled, and we see that 
the simple equations of ionization balance given in Eq. (117) need modification for actual calculations 
of the charge state distribution. This does not affect the definition of ^ given there, however. 

For iron, autoionization rates ^ lO^^-lO^"* s~^ are typical. Competing with autoionization of the 
state A** is spontaneous radiative decay. Radiative transition rates are of the same order of magnitude 
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as autoionization rates in iron, but are typically smaller by factors of a few. The case of highest prob- 
ability involves a radiative transition that fiUs the K-sheU hole. The hne energy ck is characteristic of 
the atomic structure of the ion. We write the reaction as 

A** ^ A* + eK (spontaneous radiative decay). (125) 

The product ion A* is indicated as being an excited level, although there are cases when it is a ground 
level. If the radiative transition leaves an ion in a state such that its energy lies below the first ionization 
potential, then that transition has led to radiative stabilization of the ion, i.e., the ion is no longer 
subject to autoionization. Parity-changing transitions of the type np-ls are preferred over nl-ls, where 
I ^ p. Transitions filling the K-shell hole, if n = 2, are referred to as Ka transitions, and if n = 3, K/3 
transitions. To make these concepts more concrete, we provide the following examples. 

Example 1: 11 -electron ion - autoionization 
Is22s22p63s + g ^ Is2s22p63s + g" 

(126) 

Is2s22/3s ls^2s^2p^ + e" 

The first reaction generates a photelectron, which heats the ambient plasma. The energy distribution of 
photoelectrons is continuous, and depends on the shape of the ionizing spectrum. The second reaction 
- autoionization - leaves the F-like ion in the ground level, and ejects an Auger electron. 

Example 2: 11 -electron ion - radiative stabilization by Ka emission 
ls^2s^2p^3s + e ls2s^2p^3s + e" 

Is2s22/3s ^ ls'^2s'^2p^3s -\- Ka (127) 

ls'^2s'^2p^3s Is22s22/ + cl 

The second step is a radiative stabilization of the autoionizing level, which produces a Ka photon and 

leaves an excited Ne-like ion. The latter then decays by emission of a 3s ^ 2p soft X-ray line. The end 
products are a Ne-like ground level, an Auger electron, one Ka line, and one Ne-like 3.s 2p line. 

Natural Une widths are w 0.4 A14 eV, where A denotes the sum of radiative and autoionization 
transition rates from the autoionizing level of interest, and is given as a multiple of 10^"' s^^. For a 
line with energy e, the thermal Doppler width is « 0.4 (e/6.5 keV)(/eT/100 eV)^^^ eV. Note that the 
low temperature is characteristic of X-ray photoionized plasmas. The commensurabiUty of the natural 
line width and the thermal Doppler width for iron K fluorescent lines implies large Voigt parameters 
(see, for example, [73]). Estimates of line-center optical depths based on a pure Doppler profile will 
tend to provide an overestimate by factors of a few, as well as underestimate the contribution of the 
line opacity in the damping wings. The curve of growth is affected rather dramatically; the large Voigt 
parameter results in the virtual elimination of the saturation part of the curve, leaving a transition from 
the Unear part to the damping part. 

7.2.7. KP Emission 

While not as diagnostically useful as Ka emission, the K/3 complex has the potential to provide 
coroUary information. 

Example 3: 13-electron ion - Kji emission 
\s'^2sHp^3sHp -I- e ^ \s2sHp^3sHp -|- e" 

(128) 

\s2s^2pHsHp ls^2s'^2p^3s^ + XP 
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Competing with K/3 emission in the second step are autoionization and Ka emission. There is a rela- 
tively small, but non-negligible, probabiUty that a K(3 photon is produced. Detailed calculations show 
that, for the near-neutral iron charge states, the K/3/Ka intensity ratio varies from about 0.12 (Fe II) to 
0.15 (Fe IX) [190]. 

For most cases of interest, photoionization-driven K/3 is important only when the pre-ionization 
charge state has an occupied 3p subshell in its ground configuration. To see this, we take as a coun- 
terexample, a N-like ion, with ground configuration ls^2s^2p^. While, in principle, it is possible to 
create, say, C-like ls2s^2p^3p through K-shell photoionization from the N-like ground configuration, 
this mechanism is of interest only if there is a substantial population of the excited ls^2s^2p^3p con- 
figuration, which is extremely unlikely, since the radiative lifetime of such an excited level is usually 
less than a nanosecond. Note also that it is not sufficient to have a populated 3s subshell in the ground 
configuration of the pre-ionization charge state. Decay by K/3 3s Is, because of the parity selection 
rule, is extremely improbable, since it competes with a much faster 2p ^ Is transition, as indicated in 
Example 2 above. Therefore, as a rule of thumb, K/3 is emitted only by the iron ions Fe II - Fe XIV 
(singly-ionized through Al-like) if photoionization dominates the excitation. 

One possible mechanism leading to K/3 emission in charge states more ionized then Fe XIV in- 
volves direct photoexcitation by the same radiation field responsible for creating the charge state distri- 
bution [191]. Thus for a continuum photon of energy e reactions such as ls^2s^2p^ + e ^ ls2s'^2p'^3p, 
a resonance absorption, followed by re-emission, can be efficient drivers of K/3. Similarly, again us- 
ing a C-like ion as an example, the entire Rydberg series of lines ls2s^2p^np ls^2s^2p^ can be 
energized by resonance absorption of the continuum, so that K7, Kd, etc., may appear in a spectrum. 
Gauging the overall importance of such a mechanism involves radiation transport, thereby introducing 
macroscopic properties of the gas into the calculation, which is case specific. Generally, it is found 
that the efficacy with which photoexcitation competes with photoionization in driving line emission 
decreases with plasma column density [192]. However, velocity gradients can complicate the analysis 
[193]. 

7.2.2. Fluorescence Yield 

The two-step process described above — inner-shell photoionization followed by radiative stabi- 
lization accompanied by emission of a K photon — is called fluorescence, possibly a misnomer, since 
the term applies also to other radiation processes. The Kfluorescence yield, which we denote by Yk, is 
the quotient of the rate at which K lines are generated from an irradiated sample and the rate at which 
K-shell holes are produced in the sample [194]. To motivate the concept of fluorescence yield, imag- 
ine a somewhat idealized laboratory experiment in which a beam of ions of element A interacts with a 
high-energy electron beam, say, in a crossed-beam setup. The electrons ionize the ions, increasing their 
charge by one step, and the beams are magnetically separated downstream of the interaction region. 
An experimenter counts the ionizations per unit time. Suppose also that an X-ray detector records the 
X-ray Une production, and the ratio of X-ray count rate to ionization rate is recorded. A second run 
with the beam current increased by a factor of two shows that both the ionizations per second and the 
X-ray counts per second double. Next, an ion beam of a different element B is tested with the same 
setup. It is again found that the X-ray flux doubles when the ionization rate doubles. However, the 
ratio of X rays to ionizations is different for B. One draws the conclusion that the X-ray production is 
linearly proportional to the ionization rate, but that the constant of proportionality differs from element 
to element, and, possibly, from charge state to charge state for a given element. One infers that the 
proportionality constant F is a characteristic of a given species and, specializing to K line emission, 
defines 



Yk = probability of K line emission per K shell ionization. 



(129) 
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The usual convention is to associate a fluorescence yield with the pre-ionization charge state. For ex- 
ample, Fe II K lines are associated with an Fe I yield. This makes good sense, since there are cases 
where K lines from a given charge state arise from a photoionization of an ion several charge states re- 
moved. For example, Fe III line production may add to the Fe I fluorescence yield. We also distinguish 
ennission from H-Uke and He-hke ions from that of the lower charge states. In X-ray photoionized 
plasmas, emission from H- and He-like ions results predominantly from radiative recombination into 
excited levels. One could assign a fluorescence yield to a H-like ion by taking the ratio of H-like K 
lines produced per photoionization of the H-like ion, but that is actually somewhat at odds with the 
convention used for the lower charge states. In this paper, we assume that the fluorescence yield is 
defined for neutral atoms through Li-like ions. 

In what foUows, we show that the imphcation of Eq. (129) — that F is an intrinsic atomic property 
— is false. For a set of energy levels u that lie above the ionization threshold, and a set of stabilized 
levels £, all of which belong to charge state i (see Fig. 18), one starting point for deriving an expUcit 
expression for the fluorescence yield is the emissivity, summed over all K fines of cfiarge state i, 

u t 

wfiere we assume that transitions u ^ i arc members of the Ka complex in i. Level population 
densities of level u in charge state i are denoted n,,„. Radiative transition rates are denoted by A^. The 
summed K line emissivity has c.g.s. units cm^'^ s^^. 

The level population density rii^u is found by dividing the population flux into level u by the total 
decay rate of the level, 

„ _ Xlfc '^i-i,fe/5*;M (131) 



E/lr I Aa ' 

j uj ' Z-^m um 



where is an autoionization rate connecting level u in i to level m in i -I- 1. Energy levels k in the 
pre-ionization charge state i — 1 are represented as n^-i fe. Photoionization connecting level kiai — 1 
to autoionizing levels u in i are denoted by and are calculated according to Eq. (118). 
Substituting Eq. (131) into Eq. (130) gives 

iif = ^^^rij_i,fc/3fe„r„£, (132) 

k u e 

where we define the line fluorescence yield V^i for each transition u — » £ by 



■ ui 



j uj ' L^m u 



(133) 



um 



The fine fluorescence yield depends on purely atomic quantities, and can be thought of as a radiative 
branching ratio where additional sinks are included, viz., autoionization. 

Now we "normalize" the emissivity by dividing it by the total photoionization rate, which gives an 
expression for the fluorescence yield, 

\r X]fe 

= a • (1^4) 

luk l^u ni-i,kPku 

The yield is thus a weighted average of the fine fluorescence yields, where the weightings are de- 
termined by the level population distribution of charge state i — 1 and the level-to-level photoionization 
rates, and should, in principle, be evaluated on a case-by-case basis, taking into account processes that 
affect the local excitation conditions. While complicated enough, note that Yk does not depend on the 
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level population distributions of charge states i and i + nor does it depend on the charge state distri- 
bution. We have added one further simplification, which is to ignore possible coUisional redistribution 
of levels u, which is likely to be a valid assumption for most astrophysical densities. 

The above derivation shows that a fluorescence yield cannot be considered as an intrinsic attribute 
of an atomic species, a fact that has been known for some time [195]. Equation (134) shows that Yk 
depends on the level population distribution of i — 1, which depends in turn on the electron density, 
the electron temperature, and the local radiation field. Nevertheless, tables of yields exist, and are are 
widely used in astrophysics [189] [196]. Part of the reason for this is that atomic calculations covering 
the wide range of conditions expected in astrophysical plasmas are scarce. Although this situation will 
change, in the meantime data based upon laboratory experiments, where fluorescence lines are used 
for the purposes of calibration, and "theoretical" data based upon approximation schemes, such as Z- 
scaling, continue to be used. Notwithstanding one's possible expectations that yields vary dramatically, 
depending on the plasma conditions, detailed calculations of iron fluorescence [185] [197] show that 
numerical values of Yk do not show substantial sensitivity to electron density between the low-density 
limit and the Maxwell-Boltzmann limit for the ions Fe II - Fe XVII. It is found that Yk hovers around 
the "canonical" value of 0.34 for these ions. A variation is observed for higher charge states, however. 

Since electric dipole radiative rates scale as Z"^ for An > transitions, and since Auger rates are 
roughly constant with Z, Yk increases rapidly with Z [198]. The trend of Yk with Z, based upon 
both experimental data and calculations is shown in [199]. Among the elements currently relevant to 
cosmic X-ray spectroscopy, only nickel has a higher atomic number than iron, and the nickel abundance 
is at least an order of magnitude smaller than that of iron [200]. The relatively large iron fluorescence 
yields, coupled with its high abundance, conspire to make iron K fluorescence lines the most commonly 
observed in astrophysics. Calculations of disk reflection models have shown that fluorescence lines 
from the remaining cosmically abundant elements are relatively weak [201], and this is validated by 
the lack of detection in black hole accretion disk systems, although some (e.g., siUcon Ka) are observed 
in other source classes (e.g., [202]). 

7.2.3. Energy Distribution of Iron K Lines 

To date, spectrometers used to observe extrasolar X-ray sources have, at best, resolved the iron Ka 
complex into three features: the Fe XXVI Lya doublet at 6.97 keV, the Fe XXV blend, near 6.7 keV, 
and a blend of emission from a composite of lines from Fe II-Fe XVIII near 6.4 keV. For convenience, 
we refer to the 6.4 keV blend as arising in "near-neutral" material. The energy distribution of the iron 
K Unes is shown in Fig. (19), from [185], where the potential for blending of the near-neutral Unes is 
evident. Comparing adjacent charge states, the small energy spacings of lines from near-neutral iron 
arises from the fact that changes in the screening of the atomic potential experienced by 2p electrons 
are small compared to the Ka Une energies, i.e., the potential in the n = 2 shell is dominated by 
the nuclear charge to such an extent that screening by n = 3 electrons is of small consequence. The 
blending persists into the L shell. For example, in assessing the 5 eV energy separation of the Ka line 
centroids^ of Ne-like Fe XVII and F-like Fe XVni, note that the upper configurations are ls2s^2p^3s 
and ls2s^2p®, respectively. The small difference in centroid energies results from the small screening 
of the nuclear potential by the 3s electron in the Fe XVII Is-hole state. Starting with O-like Fe XIX, 
the ion-to-ion Ka centroid separations begin to increase. The removal of a 2p electron has a relatively 
large effect on the differential screening. For example, the Ka energy centroids for Fe XVIII and Fe 
XIX are separated by 34 eV. 

The most commonly observed fluorescence Unes are those from near-neutral ions. Because of the 
blending and the small variation of Une centroid energies, the entire complex is often referred to as 



The centroid is defined as S^e^ ji /Siii, wliere e is a line energy, and j is an emissivity. 
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Fig. 19. Energy centroids for iron Ka (lower) and K/3 (upper.) Ordinate denotes thie number of bound electrons 
in the emitting ion: '2' for Fe II, etc. Figure from Kallman et al. [185]. 



"the iron line," with the understanding the "line" is likely to be a composite structure. Reference to 
table values of line energies often contain separate entries for Kai and Ka2, which partially accounts 
for the intrinsic spectral structure. The distinction between these two features results from the relative 
probabilities of a Is vacancy being filled by a radiative transition from an electron in the 2^3/2 subshell 
(Kai) or the 2pi/2 subshell {Ka2)- Thus E{Kai) > E{Ka2)- Since there are four 2^3/2 electrons 
in a filled L shell, and two 2pi/2 electrons, the Kai yield is generally given as twice the Ka2 yield. 
However, even this greatly oversimplifies the K spectra from iron ions. In fact, the Kai-Ka2 labels 
are primarily an observational convenience (when they are spectrally resolved). When one accounts 
for the myriad possibilities introduced by angular momentum coupling, and the variety of possible 
excitation conditions encountered in astrophysical plasmas, the true spectrum may consist of hundreds 
of individual lines. In some cases the Kai and Ka2 line complexes overlap in energy. 



7.2.4. Resonant Auger Destruction 

The blending of Ka Unes from near-neutral iron ions could be considered as both a blessing and a 
curse. On the one hand, the blending allows a simple approach to spectral fitting: simply add a line at or 
near 6.4 keV to the trial spectral model; no need to woiTy about a complex distribution of lines; no need 
to worry about the charge state distribution. On the other hand, blending means that while conditions 
in the disk plasma may vary considerably through the emitting region, for a large range of physical 
parameter space,, these variations are not conveyed by the spectrum. In this latter context, then, one 
might look to the L-shell ions for potential diagnostics. The Ka spectra of L-shell ions provide unique 
diagnostic information on plasmas that exist over intermediate ranges of ionization, as can be seen in 
Fig. 17. Interestingly, it appears that Ka lines from iron L-shell ions are not required in spectral fits to 
black hole accretion disk spectra, even though there is ample evidence for emission from the charge 
states that bracket them. The absence of L-shell emission has been attributed to a process known as 
resonant Auger destruction [203] [204], a radiation transfer effect that leads to severe attenuation of K 
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emission lines from L-shell ions. 

To see how resonant Auger destruction (RAD) selectively attenuates K emission from L-sheU ions, 
as opposed to M-shell ions, consider the following example. 

Example 4: 10-electron ion - resonant Auger destruction ofKa 
ls'^2s'^2p^ + e ^ Is2s22p6 + q- (photoionization) 
Is2s22p6 ^ Is22s22p5 + (Ka emission) 

(135) 

Ka + ls^2s^2p^ ls2s^2p^ (resonant absorption) 

ls2s^2p^ — > ls^2s^2p^ + e~ (Auger decay/photon destruction) 

The third step is the inverse of the second, where it is understood that the Ka photon, somewhere 
along its line of flight, encounters an ion in the same state as the product ion that resulted from the 
original Ka transition. The fourth step, autoionization, is, on average, more probable than radiative 
decay. Therefore, there is a high probability per scattering event that the photon will be destroyed. If 
the autoionizing ion in the fourth step decays radiatively instead, then, depending on the line optical 
depth of the medium, another absorption can occur, with an equally large destruction probability. The 
probability of the Ka photon surviving more than just a few scatters is very small, and very few photons 
can escape the medium. 

To make a distinction between L-shell ions and M-shell ions and their response to the RAD process, 
consider one more example. 

Example 4: 11-electron ion - an improbable resonant Auger destruction ofKa 
ls'^2s'^2p^is + Is2s22p63s + g" (photoionization) 

ls2s^2p^Zs ls^2s^2p^^s + Ka (Ka emission) 

(136) 

Ka + ls^2s22p^3s ls2s^2p^3s (resonant absorption) 

ls2s^2p^3s — > ls^2s^2p^ -|- e~ (Auger decay/photon destruction) 

While this set of reactions is possible, it is generally unimportant, since the third step is extremely 
unlikely. For this step to operate effectively, the medium would require a high optical depth in the 
relevant excited level of the ls^2s^2p^3s configuration. For the densities of interest, this will never be 
the case. By far, the dominant configuration of 10-electron ion is ls^2s^2p^, and excitation channels 
to the 2p subshell are closed. This is also the case for the RAD sequence that begins with any ion that 
has n = 3 electrons in its gound configuration. Therefore, RAD has a negUgible effect on M-sheU ions, 
but a potentially major effect on L-shell ions. 

Consider a semi-infinite slab that is irradiated from above by a hard X-ray continuum source. Only 
an upper layer of the slab corresponding roughly to unity line optical depth will contribute to Ka flu- 
orescence from L-shell ions. Below that depth, escape is prohibited by resonant Auger destruction. 
Compare this to the case of fluorescence from M-shell ions, where the fluorescing region for near- 
neutrals extends for roughly one continuum optical depth. This would suggest that the emergent Ka 
line flux from L-shell ions is dwarfed by that from the near-neutrals. This provides a plausible expla- 
nation for the absence of an L-shell component in iron Ka spectra. Moreover, in the context of spectral 
modeling of accretion disks, it appears to provide justification for simply zeroing out the Ka emission 
for L-shell ions. 



© NRC Canada 



52 



unknown Vol. 



> 0.4 




16 17 18 19 
log ion column density 



20 



0.6 



0.4 



0.2 



0.0 




16 17 18 19 20 
log ion column density 




16 17 18 19 
log ion column density 



20 




16 17 18 19 
log ion column density 



20 



Fig. 20. Effective fluorescence yields as modified by resonant Auger destruction plotted against ionic column 
density for Fe XX - Fe XXIII. Each plot shows two cases: (dotted line) level populations of low-lying states 
driven by an 80 eV half-diluted blackbody and {dashed line) no perturbing radiation field. The fluorescence yield 
of N-like Fe XX is determined by the production and escape probability of C-like Fe XXI lines, etc. 



This argument needs some modification, however. When we couch the argument only in terms 
of absorption cross-sections and ion column densities, we ignore the fact that the line optical depth 
depends on the level population of the lower level of the relevant transition, which need not be the 
ground state. For a given transition between an upper level u and a lower level I, the line optical depth, 
in terms of the hydrogen column density N^i, the elemental abundance Az, and the ionic fraction i^ion^ 
can be written 

2 

7re 

T£„ = N]iAzFiouPi fiu4>{v) (137) 

mc 

where pi is the fractional population density of level I, fiu is the absorption oscillator strength, and 
(l){v) is the line profile function in the rest frame of the absorber. Therefore, we must solve for the level 
populations for the lower levels of each Ka transition. In fact, many Ka transitions terminate on 
excited states, and such lines may be optically thin under certain conditions relevant to accretion disk 
atmospheres [205]. 

Given that the level population distribution in a charge state i partly determines the optical depth of 
Ka hnes from i, self-consistency demands that the set pi in charge state i — 1 also be calculated. The 
level populations of i — 1 determine the relative rates at which the upper levels of the Ka transitions in 
i are populated by K-shell photoionization [185]. Referring back to the discussion in §7.2.2 concerning 
the effects of local excitation conditions on Yk, this constitutes a second reason why assessments 
of RAD cannot be decoupled from microscopic considerations. It is shown in [205] that the level 
population distribution in i — 1, responding to an underlying UV/soft X-ray continuum, can significantly 
increase Yk for N-like, C-like, B-like, and Be-like iron ions. If the definiton of Yk is modified so as 
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to account for the escape probability of Ka lines from a medium of specified column density in charge 
state i, then a plot of vs. iVj allows us to evaluate the overall effect of RAD, accounting for both 
microscopic and macroscopic influences. This is illustrated in Fig. (20), which shows that RAD is not 
decisively effective in quenching Ka emission from L-shell ions until ionic column densities ^ 10^® 
cm~^ are reached. Calculations of vertical structure in irradiated accretion disks suggest that, for these 
ions, this value is most likely on the high side of what is expected in typical disk atmospheres [206]. 

7.3. X-ray Reflection 

The simplest assumption concerning vertical structure is to invoke a constant density, i.e., to forego 
a detailed consideration of the effect of the gravitational field entirely. Calculations of this type, while 
seemingly neglecting important physics, have nevertheless proven to be quite powerful in elucidating 
the physics of spectral formation in accretion disks [207] [208] [209]. With a choice for the inci- 
dent flux, the local density dictates the ionization parameter. Transport of the incident radiation down 
through the atmosphere leads to its attenuation, which results in a stratification of ^ (d^/dz is posi- 
tive). Photoionization codes are then used to generate local spectra, which, accounting for opacity in 
the overlying layers, are propagated through the atmosphere to the computational boundary, the result 
being the spectral distribution of the radiation field at the "surface." If desired, the effects of disk rota- 
tion, disk incUnation, and relativistic effects are applied, thus generating a model spectrum as observed 
at infinity. 

The first applications of constant-density models were focused on an examination of the "reflec- 
tion" of the incident continuum from cold matter, "cold" in this context meaning that H and He were 
assumed to be fuUy stripped, while the remaining elements retained all of their electrons. A fraction 
of the radiation impinging on any gas will be re-radiated, or reflected, back into the general direction 
of the radiation source. This fraction is known as the albedo. More generally, one is interested in the 
energy dependence of the albedo. High-energy photons interacting with the cold material described 
above will either be absorbed via photoionization or Compton scattered. 

A photoelectric cross-section falls roughly as above the photoionization threshold, or edge. 
Neutral iron has as its ground configuration [Ar] 3d^4s^, where the notation [Ar] symbolizes the con- 
figuration set of ground state Ar. With electrons occupying four shells, Fe I has associated with it an 
N (n = 4) edge, an M edge (n = 3), an L edge (n = 2), and a K edge (n = 1). The rapid falloff 
in the cross-sections results in the K edge dominating the total cross-section for energies above the 
7.1 keV K edge. Competing with iron opacity is Compton scattering, with a small contribution from 
photoabsorption by the remaining elements. Although the iron K edge cross-section at threshold is or- 
ders of magnitude larger than the Compton cross-section, the fact that iron is a trace element (the solar 
abundance relative to hydrogen is 4.7 x 10~^ [200]) leads to a near equality in the opacities, so that 
Compton scattering is approximately as important as the iron photoelectric opacity at the iron K edge. 
A few keV above the iron K edge, Compton scattering becomes the dominant opacity source. 

Since the typical energy lost by a photon in Compton scattering off of cold electrons is « 
E"^ jmc?, not only do high-energy photons lose more energy than low-energy photons, but they lose a 
larger fractional energy /S.E/E, as well. Compton scattering from a cold slab thus leads to a degrada- 
tion in energy of the incident spectrum, high energy electrons migrating to lower energies, with a trend 
toward a steepening relative to the incident spectral shape. At the same time, however, a fraction of the 
photons near and below ^ 10 keV are absorbed, depleting the incident flux of soft X rays. The com- 
bined effect of these two processes produces the Compton bump, an apparent excess above an incident 
power-law spectrum in the approximate range 20-50 keV [18] [19] (see Fig. 21). Accompanying this 
signature of reflection are fluorescence lines [83] [210| (see below). The overall spectrum is a sum of 
the incident power-law and the reflection spectrum. The relative contribution of the reflection spectrum 
depends on the solid angle AJ7 that the reflector subtends at the hard X-ray source [211] [212]. For 
example, a disk geometry with illumination from above gives a covering fraction AO/47r of about 1/2. 



© NRC Canada 



54 



unknown Vol. 




E (keV) 

Fig. 21. Reflection spectra for constant-density slabs for three values of ^ (as labeled), assuming normal illumi- 



nation by a power law spectrum with photon index 2. Electron density is set to 10 
set to the solar value. From Ross and Fabian [179]. 
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Some of the first serious attempts to model the vertical structure of hydrostatic accretion discs, 
including energy transport, opacity, and the equation of state, were applied to disks in cataclysmic 
variables [213] [214] [215]. These calculations made no provision for the influence of a hard X-ray 
source, however, and were not directly applicable to accreting neutron stars and black holes. Later 
developments saw the introduction of hydrostatic disk models for which the outer disk was assumed to 
be X-ray photoionized by a central source of illumination, appropriate to accreting neutron star X-ray 
binaries [216] [217] [218] [219]. 

A typical approach involves dividing the disk into a set of annuli, say, a few tens, then using a disk 
photoionization code to calculate the atmospheric structure of each annulus by integrating the equa- 
tions of hydrostatic balance and 1-D radiation transfer for a slab geometry. A common alternative is 
to treat a single annulus. A hard X-ray radiation field, incident from above at some specified angle, 
is assumed for each radius. The spectral shape is usually taken directly from observations, rather than 
calculated self-consistently. A typical choice for the relative normalizations of the X-ray field at each 
annulus is to scale it by a r^^ power law. This is motivated by the picture described earlier wherein the 
hard X rays are supposed to originate from Compton scattering of the UV/soft X-ray disk field, which 
scales as the local dissipation rate in the optically thick disk (see Eq. 8). The overall normalization 
of the X-ray flux to the UV flux is an input parameter, and is typically of order unity. In particular, 
if a single annulus is modeled, this ratio is left to the discretion of the modeler, with the assumption 
that the disk is irradiated by a localized disk flare [206] [220]. The atmosphere is divided into a set 
of zones of (possibly) variable geometrical thickness. The absorption of the incident radiation field 
leads to heating through photoelectron thermalization and Compton heating. Photoelectron mean free 
paths are sufficiently small that the energy is deposited "on the spot." A further constraint is that local 
thermal equilibrium must be satisfied, i.e., A(P, p,Fy) = 0, where the energy function depends on the 
gas pressure P, the density p and the local radiation field Fi,. The difference between heating and cool- 
ing A includes Compton scattering, bremsstrahlung cooling, photoionization heating. Auger heating. 
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collisional line cooling, and recombination cooling from H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, 
and Ni ions. Photoionization produces spectral by-products, as well, such as fluorescence lines, radia- 
tive recombination continua, and lines produced in recombination cascades. This reprocessed radiation 
contributes to the "diffuse" flux, which is added to the modified incident field, and is propagated up- 
wards and downwards to the adjacent computational ceUs. Boundary conditions at the top and bottom 
of the annulus are also prescribed. One choice [79] is to set the density at the bottom to the value given 
by the a-disk model, and to set the gas temperature at the top to the Compton temperature [73], given 
by the following integral over the radiation spectrum: 



where f/^ is the specific radiation energy density and U is the total radiation energy density. The Comp- 
ton temperature for most cases of interest is ~ 10''-10* K. 

For modeling X-ray line emission from the disk atmosphere, the commonly used assumptions of 
LTE and the diffusion approximation may not be valid. The assumption of constant density may be 
invalid, as well, since the hydrostatic equilibration time scale is small or comparable to other relevant 
timescales (cf., [69]). For calculating the discrete emission fine content of a disk spectral model, partic- 
ularly recombination emission, the vertical ionization structure and the density stratification are crucial. 
By contrast, the output iron fluorescence energy centroid varies little until the ionization reaches the L 
shell (see Fig. 18). Still, even fluorescence emission can be affected by a Compton thick, fully ionized 
gas overlying the emission sites deeper in the atmosphere [222]. Also, over certain ranges of ^, the gas 
is subject to a thermal instability that can lead to a steep transition between cold near-neutral gas and 
hot highly-ionized gas [222] [223]. 

Examples of output from a set of reflection slab models are shown in Fig. (21) [179]. In these exam- 
ples, the gas density is fixed at a constant value, and the slab is irradiated from above by a continuum 
flux with spectral shape AE~^ exp(— £J/300 keV). The value of ^ is varied by varying the normaliza- 
tion A. As described in [179], the reflection spectrum for the largest value of ^ (10'*) is nearly line-free, 
since the gas is very highly ionized. Compton broadened emission from Fe XXXV and Fe XXVI make 
the only conspicuous discrete contribution to the spectrum. At the two lower values of ^ the increas- 
ingly complex line spectra, with contributions from n = 2 ^ n = \ transitions of H-like and He-like 
ions of C, N, O, Ne, Mg, Si, S, and a few n = 3^n = 2Fe L-shell lines, are attributed to the 
progressively lower level of ionization. 

7.3.1. Thermal Instability 

X-ray irradiated gas is subject to a thermal instabihty in the 10^-10^ K temperature range [224] 
[225], suppressing X-ray line emission in that regime. Application of the Field stability criterion [226] 
indicates that a photoionized gas may become unstable when recombination cooling of H-like and He- 
hke ions is important. The T and ^ ranges where the instabihty occurs depend on the metal abundances 
and the local radiation spectrum [227], and, from the theoretical point of view, depends on the atomic 
rates used [228]. Within a range of pressure ionization parameters [175] !B = Prad/-Pgas> where Prad 
is the radiation pressure, thermal equihbrium is achieved by three distinct temperatures, only two of 
which are stable to perturbations in temperature. Physically, the instabihty causes a steep temperature 
gradient as the gas is forced to "move" between thermally stable regimes, requiring the formation 
of a transition region joining the hot and cold phases whose size could be determined by thermal 
conduction. The method by which the thermal instability is treated can significantly affect the predicted 
spectrum [79] . One issue relevant to iron K fluorescence emission is that the steep transition region may 
comprise the iron L-shell ions. Therefore, even if the effects of RAD were nil, iron L-sheU ions may 
have insufficient optical depth to contribute substantially to the overaU iron K spectrum. 
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7.5.2. Alternative X-ray Signatures of Relativistic Effects 

Prior to the launch of Chandra and XMM-Newton, the soft X-ray spectrum 1 keV) was gener- 
ally modeled as a power-law attenuated by an ionized absorber [230] - the warm absorber model. A 
radical departure from this model was suggested in [231], where it was proposed that spectral features 
observed in the 5-35 A band of the Seyfert 1 galaxies MCG-6-30-15 and Mrk 766 could be explained 
as relativistically broadened emission hues from H-like carbon, nitrogen, and oxygen. Fitting parame- 
ters - disk inclination, disk emissivity profile - were found to be roughly consistent with those found 
for the iron Ka line [27]. An analysis of a 130 ks grating observation with XMM-Newton of Mrk 766 
[232] led to the conclusion that a relativistic fine interpretation provides a better fit to the data than 
the warm absorber model. A similar conclusion was found for the Seyfert 1 galaxy NGC 4051 [233]. 
If correct, these spectral features, formed by radiative recombinaton onto bare nuclei, could provide 
additional "handles" that would allow tighter constraints to be placed on models of the accretion disk 
structure in the relativistic regime. Since the sites of formation H-like lines are spatially distinct from 
the formation sites of fluorescence lines a broader range of parameter space would be accessible. 

Still, numerous absorption fines from ionized species in the Chandra grating spectrum of MCG- 
6-30-15, also noted in [231], indicated at least some influence from a warm absorber. In fact, the 
relativistic disk interpretation in its entirety was disputed in [234]. In the model proposed there, ab- 
sorption in addition to two warm absorber components is provided by a dust component containing 
Fe02, which is embedded in the warm absorber. The warm absorber model found additional support 
from theoretical disk models, which showed that the line equivalent widths required by the relativistic 
disk interpretation far exceeded theoretical predictions [221]. An additional problem is the high nitro- 
gen to oxygen ratio required to account for the N VII feature, although a Bowen-Uke mechanism has 
been suggested as an explanation [235]. A further criticism of the relativistic disk line model is that 
reflection models predict emission from iron L-shell ions [221] that is not observed in the data. Radi- 
ation transport calculations using a Monte Carlo approach indicate no strong iron L n — > 2 emission, 
however [178]. 

In subsequent work [236] it was argued that the warm absorber model used in the Chandra inter- 
pretation of MCG-6-30-15 does not adequately describe the XMM Reflection Grating Spectrometer 
data, which has a better response at longer wavelengths. Further analysis [237] of the XMM data of 
MCG-6-30-15 showed that fits obtained with a composite model, involving both relativistic disk fines 
and a warm absorber, yield ambiguous results; if a warm absorber is added to the best-fitting disk fine 
model, then relatively little absorption is required; if relativistic lines are added to the best-fitting warm 
absorber model, then only weak emission lines are required. Given the potential importance of rela- 
tivistic recombination fine emission, there is hope that these issues will be resolved decisively in the 
near future. 

Another class of relativistic signatures has been proposed recently. The magnetic flare model [239] 
of the hot corona and hard X-ray continuum observed in accretion disks posits the existence of many 
small-scale, hot (T > 10^ K) magnetic loops. Following loop emergence, Compton cooling of loop 
material, through the interaction with softer disk photons, results in hard X-ray flaring events. Disk 
irradiation by these relatively localized flares produce fluorescence emission that is concentrated in the 
region of the disk near the flare [238]. In other words, regions of the disk that are small in both radial 
and azimuthal extent, and which radiate intense transient fluorescence emission are predicted by this 
model. The transiently fluorescing region will continue to orbit the black hole, and will emit iron fines 
that are narrow compared to the case of full disk irradiation. The centroid energy is determined by the 
gravitational redshift associated with a small range of annufi, and by the Doppler effect appropriate 
to the instantaneous projected velocities over a narrow range of azimuth. If the lifetime of such a 
region is sufficiently long, the centroid energy could be observed to change over the time scale of an 
observation. It has been shown how the fine flux received from such a "hot spot" changes with position 
along the orbit, owing to the effects of relativistic beaming and light bending [240]. Narrow spectral 
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features lying at energies below 6.4 keV have indeed been reported, e.g., [162] [241] [242]. As pointed 
out in [243], however, inferring disk parameters from these kinds of data push the Umits of current 
observational capabilities. 

8. Concluding Remarks 

The theory of accretion onto black holes has been in development for at least thirty years. Owing 
to the complex physics governing accretion, it is probably safe to predict that decades of additional 
research will be required before the theory is considered complete. Answers to outstanding questions 
concerning accretion onto black holes will not come without access to better data, and lie beyond the 
capabilities of current instrumentation. The current generation of X-ray observatories, Chandra and 
XMM-Newton are helping to identify the various avenues of inquiry [244] that will motivate future 
observations with, for example, NASA's planned Constellation-X mission and ESA's planned XEUS 
mission. 

In the meantime, on the theoretical side, referring back to Shapiro and Teukolsky's quote in §1, 
an attack on the full inhomogeneous, non-axisymmetric, time-dependent, relativistic accretion disk 
problem with coupled radiative transfer is not far beyond the horizon. In any case, it seems evident 
that disk reflection models, with their emphasis on radiation flow, and MHD simulations, with their 
emphasis on time-dependent gas flow, wiU need to be merged at some point in the not-too-distant 
future. 

Einstein developed the General Theory of Relativity at a time when it was not demanded by astro- 
nomical observations. Only Einstein's physical intuition demanded it. From Schwarzschild's solution, 
to the theory of stellar collapse, to the quasar model, to the discovery of Cyg X-1, to the supermas- 
sive black hole at the Galactic center, to the identification of relativistically modified X-ray emission 
lines, we now find that General Relativity is fully acknowledged to be an indispensable component of 
the various theoretical formalisms used to describe black hole astrophysics. One wonders if Einstein's 
skepticism regarding the existence of black holes would have remained intact in Ught of the current 
weight of evidence. 

In addition to the special and general relativity theories, Einstein's Nobel prize winning work on the 
photoelectric effect, his introduction of the A and B coefficients, his advances in statistical mechanics, 
his influence on de Broglie, Schrodinger, and Dirac, and many other of his contributions to theoretical 
physics, figure prominently in modern research on the topics discussed in this paper. 

Acknowledgements 

We thank Gordon Drake and Peter Beiersdorfer for their advice and encouragement. Our apprecia- 
tion to Chris Mauche, whose careful reading of the text led to improvements in the overall presentation. 
Thanks also to Rob Fender, Mario Jimenez-Garate, Tim Kalhnan, Ben Mathiesen, Paul Nandra, Mar- 
tin Pessah, Randy Ross, and Masao Sako for contributions. An anonymous referee contributed several 
helpful suggestions. Work at LLNL was performed under the auspices of the U.S. Department of 
Energy by the University of California Lawrence Livermore National Laboratory under contract No. 
W-7405-Eng-48. 

References 

1. A. Pais, Subtle is the Lord (Oxford University Press; 1982). 

2. R. Giacconi, H. Gursky, F.R. Paolini, and B.B. Rossi, Phys. Rev. Lett. 9, 439 (1962). 

3. M. Schmidt, Nature 197, 1040 (1963). 

4. A. Hewish, S.J. Bell, J.D.H. Pilkington, PR Scott, and R.A. Collins, Nature 217, 709 (1968). 

5. M. Nauenberg and G.J. Chapline, Ap. J. 179, 277 (1973). 

6. C.E. Rhoades and R. Ruffini, Phys. Rev. Lett. 32, 324 (1974). 



© NRC Canada 



58 



imkiiown Vol. 



7. I.D. Novikov, in Theory of Black Hole Accretion Disks, eds., M.A. Abramowicz, G. Bjomsson, and J.E. 
Pringle (Cambridge University Press), p. 41 (1998). 

8. C.W. Allen, Astrophysical Quantities, Second Edition (The Athlone Press, University of London; 1963). 

9. A.C. Fabian, W.N. Brandt, R.G. McMahon, and l.M. Hook, M.N.R.A.S. 291, L5 (1997). 

10. S.L. Shapiro and S.A. Teukolsky, Black Holes, White Dwarfs, and Neutron Stars (John Wiley and Sons; 
1983). 

11. M. Elvis, et al., M.N.R.A.S. 183, 129 (1978). 

12. M. Oda, P. Gorenstein, H. Gursky, E. Kellogg, E. Schreier, H. Tananbatim, and R. Giacconi, Ap. J. 166, 

LI (1971). 

13. M. van der Klis, in X-Ray Binaries, ed., W.H.G. Lewin, J. van Paradijs, and E.P.J, van den Heuvel 
(Cambridge University Press: Cambridge), p. 252 (1995). 

14. H. Kuneida, T.J. Turner, H. Awaki, K. Koyama, R.F. Mushotzky, and Y. Tsusaka, Nature 345, 786 
(1990). 

15. R.R Mushotzky, et al., Ap. J. 235, 377 (1980). 

16. N.E. White, Adv Space Res. 3, 9 (1984). 

17. M.M. Basko, Ap. J. 223, 268 (1978). 

18. RW. Guilbert and M.J. Rees, M.N.R.A.S. 233, 475 (1988). 

19. A.R Lightman and T.R. White, Ap. J. 233, 57 (1988). 

20. K. Nandra, in X-ray Astronomy: Stellar Endpoints, AGN, and the Diffuse X-ray Background, eds., N. 
White, G. Malaguti, and G. Palumbo (American Institute of Physics), p. 264 (2001). 

21. A.C. Fabian, M.J. Rees, L. Stella, and N.E. White, M.N.R.A.S. 238, 729 (1989). 

22. K. Chen, J.R Halpem, and A.V. Filippenko, Ap. J. 339, 742 (1989). 

23. R Barr, N.E. White, and C.G. Page, M.N.R.A.S. 216, 65 (1985). 

24. H. van der Woerd, N.E. White, and S.M. Kahn, Ap. J. 344, 320 (1989). 

25. T.R. Kallman and N.E. White, Ap. J. 341, 955 (1989). 

26. Y. Tanaka, et al.. Nature 375, 659 (1995). 

27. K. Nandra, l.M. George, R.F. Mushotzky, T.J. Turner, and T. Yaqoob, Ap. J. 477, 602 (1997). 

28. R.F Mushotzky, C. Done, and K.A. Pounds, ARA&A 31, 717 (1993). 

29. M. Elvis, Ap. J. 545, 63 (2000). 

30. L. Titarchuk, D. Kazanas, and PA. Becker, Ap. J. 598, 411 (2003). 

31. M. van der Klis, Ann. Rev. Astron. Astrophys. 38, 717 (2001). 

32. J. McClintock and R.A. Remillard, in Compact Stellar X-ray Sources, eds., W.H.G. Lewin and M. van 
der Klis (Cambridge University Press; 2003). 

33. A.C. Fabian, K. Iwasawa, C.S. Reynolds, and A.J. Young, PA.S.P 112, 1145 (2000). 

34. C.S. Reynolds and M.A. Nowak, Phys. Rep. 377, Issue 6, 389 (2003). 

35. J.R. Oppenheimer and G.M. Volkoff, Phys. Rev. 55, 374 (1939). 

36. J.R. Oppenheimer and H. Snyder, Phys. Rev. 56, 455 (1939). 

37. F Hoyle and R.A. Lyttleton, Proc. Camb. Phil. Soc. 35, 405 (1939). 

38. A. Treves, L. Maraschi, and M. Abramowicz, Accretion: A Collection of Influential Papers (World 
Scientific Pubhshing Co. Pte. Ltd; 1989). 

39. J. Frank, A. King, and D. Raine, Accretion Power in Astrophysics, Second Edition, (Cambridge: 
Cambridge University Press; 1992). 

40. O. Blaes, in Accretion Disks, Jets, and High Energy Phenomena in Astrophysics, ed., V. Beskin, et al., 
Les Houches Summer School, 78, 137 (2004). 

41. K. Davidson and J.R Ostriker, Ap. J. 179, 585 (1973). 

42. S.L. Shapiro, Ap. J. 180, 531 (1973). 

43. V.F Shvartsman, Soviet Astr. 15, 37 (1971). 

44. R Meszaros, A&A 44, 59 (1975). 

45. S.H. Lubow and F.H. Shu, Ap. J. 198, 383 (1975). 

46. F. Verbunt and E.P.J, van den Heuvel, in X-Ray Binaries, ed., W.H.G. Lewin, J. van Paradijs, and E.P.J, 
van den Heuvel (Cambridge University Press: Cambridge), p. 457 (1995). 

47. R.P Kraft, Ap. J. 135, 408 (1961). 

48. J.E. Pringle and M.J. Rees, A&A 21, 1 (1972). 

49. N.I Shakura and R.A. Sunyaev, A&A 24, 337 (1973). 



© NRC Canada 



Liedahl and Torres 



59 



50. M.A. Abramowicz, B. Czemy, J.-P. Lasota, and E. Szuszkiewicz, Ap. J., 332, 646 (1988). 

51. X. Chen and R.E. Taam, Ap. J. 412, 254 (1993). 

52. C. Luo and E.P. Liang, E.P., Ap. J. 498, 307 (1998). 

53. S. Ichimaru, Ap. J. 214, 840 (1977). 

54. R. Narayan and 1. Yi, Ap. J. 452, 710 (1995). 

55. Y. Tanaka and W.H.G. Lewin, in X-Ray Binaries, ed., W.H.G. Lewin, J. van Paradijs, and E.P.J, van den 
Heuvel (Cambridge University Press: Cambridge), p. 126 (1995). 

56. Y. Tanaka and N. Shibazaki, ARA&A, 34, 607 (1996). 

57. R.D. Blandford and M.C. Begelman, M.N.R.A.S. 303, 1 (1999). 

58. E. Quataert, E. and A. Gruzinov, Ap. J. 539, 809 (2000). 

59. M.A. Malkan, Ap. J. 268, 582 (1983). 

60. J.H. Krolik, in Theory of Black Hole Accretion Disks, ed. M.A. Abramowicz, G. Bjomsson, and J.E. 
Pringle (Cambridge University Press), p. 134 (1998). 

61. X. Chen, M.A. Abramowicz, J. Lasota, R. Narayan, and I. Yi, Ap. J., 443, L61 (1995). 

62. S.A. Balbus and J.R Hawley, Rev. Mod. Phys. 70, 1 (1998). 

63. J.F. Hawley, C.R Gammie, and S.A. Balbus, Ap. J. 440, 742 (1995). 

64. G.A. Shields, Nature 272, 706 (1978). 

65. M.A. Malkan and W.L.W. Sargent, Ap. J. 254, 22 (1982). 

66. M. Elvis, et al., Ap. J. SuppL, 95, 1 (1994). 

67. B. Czemy and M. Elvis, Ap. J. 321, 305 (1987). 

68. N.E. White, L. Stella, and A.N. Parmar, Ap. J. 324, 363 (1988). 

69. R.R. Ross, A.C. Fabian, and S. Minshige, M.N.R.A.S. 258, 189 (1992). 

70. L. Maraschi and S. Molendi, Ap. J. 353, 452 (1990). 

71. J. Poutanen, J.H. Krolik, and R Ryde, M.N.R.A.S., 292, L21 (1997). 

72. A.S. Kompaneets, Soviet Physics JETP 4, 730 (1957). 

73. G.B. Rybicki and A.P. Lightman, Radiative Processes in Astrophysics (Wiley-Interscience; 1979). 

74. P.J.E. Peebles, Principles of Physical Cosmology (Princeton University Press; 1993). 

75. G.S. Bistnovatyii-Kogan, Y.B. Zeldovich, and R.A. Sunyaev, Sov. Astron. AJ 15, 17 (1971). 

76. A.A. Zdziarski, Ap. J. 289, 514 (1985). 

77. R. Svennson, in Radiation Hydrodynamics in Stars and Compact Objects, lAU Coll. No. 89, eds. D. 
Mihalas and K.-H. Winkler (Springer- Verlag), p. 325 (1986). 

78. G. Ghisellini and R Haardt, Ap. J. 429, L53 (1994). 

79. M.A. Jimenez-Garate, J.C. Raymond, D.A. Liedahl, and C.J. Hailey, Ap. J. 558, 448 (2001). 

80. D.A. Liedahl, in X-Ray Spectroscopy in Astrophysics, eds J. van Paradijs and J. Bleeker, Springer 
Lecture Notes in Physics, p. 189 (1999). 

81. M.A. Bautista, T.R. Kallman, L. Angelini, D.A. Liedahl, and D.R Smits, Ap. J. 509, 848 (1998). 

82. J.B. Kingdon and G.J. Rerland, Ap. J. 516, L107. 

83. l.M. George and A.C. Fabian, M.N.R.A.S. 249, 352 (1991). 

84. Q.Z. Liu, J. van Paradjis, and E.PJ. van den Heuvel, A&A Suppl. 147, 25 (2000). 

85. Q.Z. Liu, J. van Paradjis, and E.PJ. van den Heuvel, A&A 368, 1021 (2000). 

86. P. Charles, in Theory of Black Hole Accretion Disks, ed. M.A. Abramowicz, G. Bjomsson, and J.E. 
Pringle (Cambridge University Press), p. 1 (1998). 

87. C.T. Bolton, Nature, 235, 271 (1972). 

88. E.P.J, van den Heuvel, in ESA, Environment Observation and Climate Modelling Through International 
Space Projects (SEE N93-2387808-88), (1992). 

89. G.E. Brown and H.A. Bethe, Ap. J. 423, 659 (1994). 

90. RX. Timmes, S.E. Woosley, and T.A. Weaver, Ap. J. 457, 834 (1996). 

91. E. Agol, M. Kamionkowski, L.V.E. Koopmans, and R.D. Blandford, Ap. J. 576, L131 (2002). 

92. M.R. Garcia, J.E. McClintock, R. Narayan, R Callanan, D. Barret, and S.S. Murray, Ap. J. 553, 47 
(2001). 

93. A.M. Stirling, et al., M.N.R.A.S., 327, 1273 (2001). 

94. R.P Render, M.N.R.A.S. 322, 377 (2001). 

95. G.E. Romero, M.M. Kaufman Bemado, and I.E. Mirabel, A&A, 393, L61 (2002). 

96. D.R Torres, G.E. Romero, J. A. Combi, P. BenagUa, H. Andemach, and B. Punsly, A&A 370, 468 



© NRC Canada 



60 



imkiiown Vol. 



(2001). 

97. D.F. Torres, G.E. Romero, X. Barcons, and Y. Lu, Ap. J., in press, astro-ph/0503186 (2005). 

98. J.B. Dove, J. Wilms, M. Maisack, and M.C. Begelman, Ap. J., 487, 759 (1997). 

99. A.A. Esin, J.E. McClintock, and R. Narayan, Ap. J., 489, 865 (1997). 

100. A.A. Esin, R. Narayan, W. Cui, J.E. Grove, and S. Zhang, Ap. J., 505, 854 (1998). 

101. R.P Fender, T.M. Belloni, and E. Gallo, M.N.R.A.S. 335, 1105 (2004). 

102. J.M. Miller, et al., Ap. J., 578, 348 (2002). 

103. A. Martocchia, et al, A&A 387, 215 (2002). 

104. J.M. Miller, et al., Ap. J. 577, L15 (2002). 

105. J.M. Miller, et al., Ap. J. 570, L69 (2002). 

106. J.H. Krolik, AGN: From the Central Black Hole to the Galactic Environment (Princeton University Press; 
1999). 

107. L. Woltjer, in Active Galactic Nuclei - Saas-Fee Advanced Course 20, eds., T.J.-L. Courvoisier and M. 
Mayor (Springer- Verlag), p. 1 (1990). 

108. R. Antonucci, Ann. Rev Astron. Astrophys. 31, 473 (1993). 

109. H. Netzer, in Active Galactic Nuclei - Saas-Fee Advanced Course 20, eds., T.J.-L. Courvoisier and M. 
Mayor (Springer- Verlag), p. 57 (1990). 

110. R.F. Mushotzky, in New Visions of the X-ray Universe, astr-pli/0203310 (2002). 

111. R.E Mushotzky, L.L. Cowie, A.J. Barger, and K.A. Amaud, Nature 404, 459 (2000). 

112. G. Hasinger, et al., A&A 365, 45 (2001). 

113. M. Miyoshi, J. Moran, J. Hermstein, L. GreenhiU, N. Nakai, P. Diamond, and M. Inoue, Nature 373, 127 
(1995). 

114. L. Ferrarese L. and H. Ford, Space Science Reviews 116, 523 (2005). 

115. J. Magorrian, et al., A.J. 115, 2285 (1998). 

116. J. Kormendy and K. Gebhardt, in 20th Texas Symposium on Relativistic Astrophysics, AIP Conference 
Proceedings, Vol. 586, ed., J.C. Wheeler, p. 363 (2001). 

117. F. Melia and H. Falcke, Ann. Rev Astron. Astrophys. 39, 309 (2001). 

118. R. Schodel, et al.. Nature 419, 694 (2002). 

119. A.M. Ghez, Ap. J. 586, L127 (2003). 

120. A.M. Ghez, S. Salim, S.D. Homstein, A. Tanner, M. Morris, E.E. Becklin, and G. Duchene, Ap. J. 620, 
744 (2005). 

121. D.F Torres, S. Capozziello, and G. Lambiase, Phys. Rev D62, 104012 (2000). 

122. J. Kormendy and D. Richstone, Ann. Rev. Astron. Astrophys. 33, 581 (1995). 

123. J. Wilms, C.S. Reynolds, M.C. Begelman, J. Reeves, S. Molendi, R. Staubert, and E. Kendziorra, 
M.N.R.A.S. 328, L27 (2001). 

124. A.C. Fabian, et al., M.N.R.A.S. 335, LI (2002). 

125. K.O. Mason, G. Branduardi-Raymont, P.M. Ogle, M.J. Page, E.M. Puchnarewicz, and N.J. Salvi, Adv. 
Space Res. 34, 2610 (2004). 

126. K. Nandra, l.M. George, R.F. Mushotzky, T.J. Turner, and T. Yaqoob, Ap. J. 523, L17 (1999). 

127. C.S. Reynolds, in Probing the Physics of Active Galactic Nuclei, ASP Conference Proceedings, Vol. 224, 
eds., B.M. Peterson, R.W Pogge, and R.S. Polidan, p. 105 (2001). 

128. I.R Mirabel, et al., Nature 358, 215 (1992). 

129. I.E. Mirabel and L.R Rodriguez, Nature 392, 673 (1998). 

130. R.P. Fender and T. Maccarone, in Cosmic Gamma-ray Sources, eds., K.S. Cheng, et al., p. 205 (2003). 

131. B.J. Sams, A. Eckart, and R. Sunyaev, Nature 392, 673 (1998). 

132. M. Rees, Nature 211, 468 (1966). 

133. I.R Mirabel, A&A 330, L9 (1998). 

134. H. van der Laan, Nature 211, 1131 (1966). 

135. A.R Marscher et al.. Nature 417, 625 (2002). 

136. CM. MiUer and E.J.M. Colbert, International Journal of Modern Physics D 13, 1 (2004). 

137. S. Weinberg, Gravitation and Cosmology (John Wiley and Sons; 1972). 

138. B.F Schutz, A First Course in General Relativity (Cambridge University Press; 1985). 

139. L. Stella, Nature 344, 747 (1990). 

140. A. J. Young, R.R. Ross, and A.C. Fabian, M.N.R.A.S 300, Lll (1998). 



© NRC Canada 



Liedahl and Torres 



61 



141. C.S. Reynolds, A.J. Young, M.C. Begelman, and A.C. Fabian, Ap. J. 514, 164 (1999). 

142. J.C. Lee, A.C. Fabian, C.S. Reynolds, W.N. Brandt, and K. Iwasawa, M.N.R.A.S. 318, 857 (2000). 

143. J. Chiang, et al., Ap. J. 528, 292 (2000). 

144. C.S. Reynolds, Ap. J. 533, 811 (2000). 

145. G. Miniutti and A.C. Fabian, M.N.R.A.S. 349, 1435 (2004). 

146. J.M. Bardeen, W.H. Press, and S.A. Teukolsky, Ap. J. 178, 347 (1972). 

147. J.M. Bardeen, Nature 226, 64 (1970). 

148. R.R Kerr, Phys. Rev. Lett. 11, 237 (1963). 

149. K.S. Thome, Black Holes and Time Warps: Einstein's Outrageous Legacy (W.W. Norton and Company; 
1994). 

150. R.H. Boyer and R.W. Lindquist, J. Math. Phys. 8, 265 (1967). 

151. K.S. Thome, Ap. J. 191, 507 (1974). 

152. E.F. Taylor and J.A. Wheeler, Exploring Black Holes: Introduction to General Relativity (Addison 
Wesley Longman; 2000). 

153. R. Penrose, Riv. Nuovo Cimento, 1, 252 (1969). 
154. 1. Ciufolini and B.C. Pavlis, Nature 431, 958 (2004). 

155. 1.D. Novikov and K.S. Thome, in Black Holes, ed., C. DeWitt and B. DeWitt (Gordon and Breach: New 
York; 1972). 

156. A.P. Lightman, W.H. Press, R.H. Price, and S.A. Teukolsky, Problem Book in Relativity and Gravitation 
(Princeton University Press; 1975). 

157. B. Aschenbach, A&A 425, 1075 (2004). 

158. A. Laor, Ap. J. 376, 90 (1991). 

159. K.A. Amaud, in Astronomical Data Analysis Software and Systems V, ASP Cof. Ser. Vol. 101, p. 17, 
eds., G.H. Jacoby and J. Barnes (1996). 

160. J.C. Houck, High Resolution X-ray Spectroscopy with XMM-Newton and Chandra, ed., G. Branduardi- 
Raymont, published electronically and stored on CD (2002). 

161. K. Iwasawa, et al., M.N.R.A.S. 282, 1038 (1996). 

162. K. Iwasawa, A.C. Fabian, A.J. Young, H. Inoue, and C. Matsumoto, M.N.R.A.S. 306, L19 (1999). 

163. Y. Dabrowski, A.C. Fabian, K. Iwasawa, A.N. Lasenby, and C.S. Reynolds, M.N.R.A.S. 288, Lll 
(1997). 

164. C.S. Reynolds and M.C. Begelman, Ap. J. 488, 109 (1997). 

165. J.H. Krolik and J.F Hawley, Ap. J. 573, 754 (2002). 

166. K. Rauch and R.D. Blandford, Ap. J. 421, 46 (1994). 

167. K. Beckwith and C. Done, M.N.R.A.S. 352, 353 (2004). 

168. M. Dovciak, V. Karas, and T. Yaqoob, Ap. J. Suppl. 153, 205 (2004). 

169. L.-X. Li, E.R. Zimmerman, R. Narayan, and J.E. McCUntock, Ap. J. Suppl. 157, 335 (2005). 

170. C.W. Misner, K.S. Thorne, and J.A. Wheeler, Gravitation (W.H. Freeman and Co., San Francisco; 1973). 

171. A.C. Fabian, et al., M.N.R.A.S. 277, Lll (1995). 

172. G.J. Ferland, K.T. Korista, D.A. Vemer, J.W. Ferguson, J.B. Kingdon, and E.M. Vemer, P.A.S.P. 110, 
761 (1998). 

173. T.R. Kallman and J.H. Krolik, XSTAR vl.43, HEASARC (NASA/GSFC), Greenbelt, MD (1999). 

174. J.H. Krolik and T.R. Kallman, Ap. J., 286, 366 (1984). 

175. J.K. Krolik, C.F McKee, and C.B. Tarter, Ap. J. 249, 422 (1981). 

176. 1. Hubeny, in Spectroscopic Challenges of Photoionized Plasmas, ASP Conference Series, Vol. 197, p. 
175, eds., G. Ferland and D. Savin (2001). 

177. A.-M. Dumont and S. Collin, in Spectroscopic Challenges of Photoionized Plasmas, ASP Conference 
Series, Vol. 197, p. 231, eds., G. Ferland and D. Savin (2001). 

178. C.W. Mauche, D.A. Liedahl, B.F. Mathiesen, M.A. Jimenez-Garate, and J.C. Raymond, Ap. J. 606, 168 
(2004). 

179. R.R. Ross and A.C. Fabian, M.N.R.A.S. 358, 211 (2005). 

180. K. Davidson and H. Netzer, Rev. Mod. Phys. 51, 715 (1979). 

181. J.R Halpem and J.E. Grindlay, Ap. J. 242, 1041 (1980). 

182. T.R. Kalhnan and R. McCray, Ap. J. Suppl. 50, 263 (1982). 

183. T.R. Kallman, in Spectroscopic Challenges of Photoionized Plasmas, ASP Conference Series, Vol. 247, 



© NRC Canada 



62 



imkiiown Vol. 



p. 175, eds., G. Ferland and D. Savin (2001). 

184. C.B. Tarter, W. Tucker, and E.E. Salpeter, Ap. J. 156, 943 (1969). 

185. T.R. Kallman, P. Palmeri, M. Bautista, C. Mendoza, and J.H. Krolik, Ap. J. Suppl. 155, 675 (2004). 

186. K.A. Nandra, K.A. Pounds, G.C. Stewart, A.C. Fabian, and M.J. Rees, M.N.R.A.S. 236, 39 (1989). 

187. M. Matsuoka, L. Piro, M. Yamauchi, and T. Murakami, Ap. J. 361, 440 (1990). 

188. J.C. Weislieit, Ap. J. 190, 735 (1974). 

189. J.S. Kaastra and R. Mewe, A&A Suppl. 97, 443 (1993). 

190. R Palmeri, C. Mendoza, T.R. Kallman, M.A. Bautista, and M. Melendez, A&A 410, 359 (2003). 

191. D.A. Liedahl, Atomic Data for X-ray Astronomy, 25th Meeting of the lAU, Joint Discussion 17, 22 July 

2003, Sydney, in press (2005). 

192. A. Kinkhabwala, et al., Ap. J. 575, 732 (2002). 

193. RS. Wojdowski, D.A. Liedahl, M. Sako, S.M. Kahn, and F. Paerels, Ap. J. 582, 959 (2003). 

194. W. Bambynek, et al.. Rev. Mod. Rhys. 44, 716 (1972). 

195. V.L. Jacobs, G.A. Doschek, J.F. Seely, and R.D. Cowan, Rhys. Rev. A 39, 2411 (1989). 

196. L.L. House, Ap. J. Suppl. 155, 21 (1969). 

197. R Palmeri, C. Mendoza, T.R. Kallman, and M.A. Bautista, A&A 403, 1175 (2003). 

198. R.D. Cowan, The Theory of Atomic Structure (Berkeley, Los Angeles, London: University of CaUfomia 
Press), p. 589 (1981). 

199. M.O. Krause, J. Phys. Chem. Ref. Data 8, 307 (1979). 

200. E. Anders and N. Grevesse, Geochimica et Cosmochimica Acta 53, 197 (1989). 

201. G. Matt, A.C. Fabian, and C.S. Reynolds, M.N.R.A.S 289, 175 (1997). 

202. M. Sako, D.A. Liedahl, S.M. Kahn, and F. Paerels, Ap. J. 525, 921 (1999). 

203. R.R. Ross, A.C. Fabian, and W.N. Brandt, M.N.R.A.S. 278, 1082 (1996). 

204. D.L. Band, R.l. Klein, J.l. Castor, and J.K. Nash, Ap. J. 362, 90 (1990). 

205. D.A. Liedahl, in X-ray Diagnostics of Astrophysical Plasmas, ed., R.K. Smith, (American Institute of 
Physics), p. 99 (2005). 

206. S. Nayakshin and T.R. and Kallman, Ap. J. 546, 406 (2001). 

207. R.R. Ross and A.C. Fabian, M.N.R.A.S. 261, 74 (1993). 

208. G. Matt, A.C. Fabian, and R.R. Ross, M.N.R.A.S. 262, 179 (1993). 

209. RT. Zycki, J.H. Krolik, A.A. Zdziarski, and T.R. Kallman, Ap. J. 437, 597 (1994). 

210. K.A. Pounds, K. Nandra, G.C. Stewart, I.M. George, and A.C. Fabian, Nature 344, 132 (1990). 

211. G. Matt, G.C. Perola, and L. Piro, A&A 245, 75 (1991). 

212. G. Matt, G.C. Perola, L. Piro, and L. Stella, A&A 257, 63 (1992). 

213. F. Meyer and E. Meyer-Hofmeister, A&A 106, 34 (1982). 

214. S. Mineshige and Y. Osaki, P.A.S.J. 35, 377 (1983). 

215. J.K. Canizzo and J.C. Wheeler, Ap. J. Suppl. 40, 1 (1984). 

216. Y.-K. Ko and T.R. Kallman, Ap. J. 374, 721 (1991). 

217. Y.-K. Ko and T.R. Kallman, Ap. J. 431, 273 (1994). 

218. J.C. Raymond, Ap. J. 412, 267 (1993). 

219. A. Rozanska and B. Czemy, Acta Astronomica 46, 233 (1996). 

220. D.R. Ballantyne and R.R. Ross, M.N.R.A.S. 332, 777 (2002). 

221. D.R. Ballantyne, R.R. Ross, and A.C. Fabian, M.N.R.A.S. 336, 867 (2002). 

222. S. Nayakshin, D. Kazanas, and T.R. Kallman, Ap. J. 537, 83 (2000). 

223. M.A. Jimenez-Garate, J.C. Raymond, and D.A. Liedahl, Ap. J. 581, 1297 (2002). 

224. J. Buff and R. McCray, Ap. J. 189, 147 (1974). 

225. G.B. Field, D.W. Goldsmith, and H.J. Habing, Ap. J. 155, L149 (1969). 

226. G.B. Field, Ap. J. 142, 531 (1965). 

227. C.J. Hess, S.M. Kahn, and RB.S. Paerels, Ap. J. 478, 94 (1997). 

228. D.W. Savin, et al., Ap. J. Suppl. 123, 687 (1999). 

229. G. Matt, A.C. Fabian, and R.R. Ross, M.N.R.A.S. 280, 823 (1996). 

230. J.R Halpern, Ap. J. 281, 90 (1984). 

231. G. Branduardi-Raymont, M. Sako, S.M. Kahn, A.C. Brinkman, J.S. Kaastra, and M.J. Page, A&A 365, 
L140 (2001). 

232. K.O. Mason, et al., Ap. J. 582, 95 (2003). 



© NRC Canada 



Liedahl and Torres 



63 



233. P.M. Ogle, K.O. Mason, M.J. Page, N.J. Salvi, F.A. Cordova, I.M. McHardy, and W.C. Priedhorsky, Ap. 
J. 606, 151 (2004). 

234. J.C. Lee. et al., Ap. J. 554, L13 (2001). 

235. M. Sako, Ap. J. 594, 1108 (2003). 

236. M. Sako, et al., Ap. J. 596, 114 (2003). 

237. A.K. Turner, A.C. Fabian, S. Vaughan, and J.C. Lee, M.N.R.A.S. 346, 833 (2003). 

238. S. Nayakshin and D. Kazanas, Ap. J. 553, L141 (2001). 

239. A.A. Galeev, R. Rosner, and G.S. Vaiana, Ap. J. 229, 318 (1979). 

240. M. Dovciak, S. Bianchi, M. Guainazzi, V. Karas, and G. Matt, M.N.R.A.S. 350, 745 (2004). 

241. T.J. Turner, et al., Ap. J. 574, L123 (2002). 

242. T. Yaqoob, I.M. George, T.R. Kallman, U. Padmanabhan, K.A. Weaver, and T.J. Turner, Ap. J. 596, 85 
(2003). 

243. K. Iwasawa, G. Miniutti, and A.C. Fabian, M.N.R.A.S. 355, 1073 (2004). 

244. FB.S. Paerels and S.M. Kahn, ARA&A 41, 291 (2003). 



© NRC Canada 



